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ABSTRACT 



Opacity effects in relativistic sources of high-energy gamma-rays, such as 
gamma-ray bursts (GRBs) or Blazars, can probe the Lorentz factor of the out- 
flow as well as the distance of the emission site from the source, and thus help 
constrain the composition of the outflow (protons, pairs, magnetic field) and the 
emission mechanism. Most previous works consider the opacity in steady state. 
Here we study the effects of the time dependence of the opacity to pair production 
(/y/y _^ e + e~) in an impulsive relativistic source, which may be relevant for the 
prompt gamma-ray emission in GRBs or flares in Blazars. We present a simple, 
yet rich, semi-analytic model for the time and energy dependence of the optical 
depth, T"ry<y - 111 which a thin spherical shell expands ultra-relativistically and emits 
isotropically in its own rest frame over a finite range of radii, Rq < R < Rq + AR. 
This is particularly relevant for GRB internal shocks. We find that in an impul- 
sive source (AR < R ), while the instantaneous spectrum (which is typically 
hard to measure due to poor photon statistics) has an exponential cutoff above 
the photon energy S\(T) where t 77 (£i) = 1, the time integrated spectrum (which 
is easier to measure) has a power-law high-energy tail above the photon energy 
£1* ~ £i(AT) where AT is the duration of the emission episode. Furthermore, 
photons with energies e > Eu are expected to arrive mainly near the onset of the 
spike in the light curve or flare, which corresponds to the short emission episode. 
This arises since in such impulsive sources it takes time to build-up the (target) 
photon field, and thus the optical depth r 77 (e) initially increases with time and 
£\(T) correspondingly decreases with time, so that photons of energy e > Eu are 
able to escape the source mainly very early on while £\(T) > e. As the source ap- 
proaches a quasi-steady state (AR ^> Rq), the time integrated spectrum develops 
an exponential cutoff, while the power-law tail becomes increasingly suppressed. 
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1. Introduction and motivation 

Astrophysical sources of gamma-rays that are both compact and very luminous may 
be optically thick to pair production (77 — > e + e~) within the source. The corresponding 
optical depth, r 77 , is usually an increasing function of the photon energy, and therefore 
a large optical depth would prevent the escape of high-energy photons from the source, 
causing a high-energy cutoff in the observed spectrum. For sufficiently high optical depths, 
enough e + e~ pairs may be produced, so that the optical depth of all photons (even low energy 
photons that are optically thin to pair production) to scattering on these electrons/positrons 
would be much larger than unity, in which case the photon energy spectrum would be 
thermalized. The size of the gamma-ray emitting region is usually hard to constrain directly 
from observations, since the angular resolution of gamma-ray telescopes is much poorer than 
their counterparts in lower energy photons (e.g. X-rays, optical, or radio). Nevertheless, the 
physical properties of the emitting region can be constrained using compactness arguments, 
and the observed properties of the source. In particular, rapid flux variability of the source 
is often used in order to set upper limits on the size of the emitting region, making highly 
variable sources with significant non-thermal high-energy emission a prime target for such 
analysis. One of the best examples for such sources are gamma-ray bursts (GRBs), and 
we shall focus on them below, although most of our analysis has a much broader range of 
applicability (similar opacity consider ations have also been used to co nstrain the properties 



of other sources, such as Blazars, e.g. ISikora. Begelman fc Reesl 1 19941 1 



It has been realized early on that, in GRBs, pair production within the s ource is expe cted 



to cause a high-energy cutoff in the observed photon energy spectrum (see lPiranll2005l . and 
references therein). Naively, if the source shows significant flux variability on an observed 
time scale of AT, its size is inferred to be R < cAT/(l + z) where z is its cosmological 
redshift, and the optical depth to pair production at a dimensionless photon energy e = 
E ph /m e c 2 is r 77 (e) ~ <7 T L 1/£ ( 1+z) /4:7im e c 3 R > a T L 1/E{1+z) (l + z)/47rm e c /L AT ~ 10 14 (1 + 
z)[L 1/e(1+z) /(10 51 ergs- 1 )][AT/(l ms)]" 1 , where L E = F e/(1+t) 4n<P L (l + z)- 1 and F e are 
the source isotropic equivalent luminosity and observed flux per unit dimensionless photon 
energy, and d^ is the luminosity distance to the source. For GRBs the (observed part of 
the) eF £ spectrum typically peaks around e ~ 1, and being at cosmological distances their 
isotropic equivalent luminosity is typically in the range of 10 50 — 10 53 erg s -1 . Furthermore, 
they often show significant variability down to millisecond timescales. This implies huge 
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values of r 77 , as high as ~ 10 , under the above naive assumptions. Such huge optical 
depths are clearly inconsistent with the non-thermal GRB spectrum, wh ich has a significa nt 
power law high-energy tail. This is known as the compactness problem (IRudermanl 1 1 9 75l ) . 



If the source is moving relativistically toward us with a Lorentz factor T ^> 1, then 
in its own rest frame the photons have smaller energies, e' ~ e(l + z)/T, while in the lab 
frame (i.e. the rest frame of the central source) most of the photons angles 
< 1/r relative to its direction of motion. The latter implies that in the lab frame the 
typical angle between the directions of the interacting photons is #12 ~ 1/r, which has two 
effects. First, it increases the threshold for pair production, (1 + z) 2 eiE2 > 2/(1 — cos #12), 
to (1 + z) 2 eiE2 > T 2 (compared to e' x e' 2 > 1 for the roughly isotropic distribution of angles 
between the directions of the interacting photons in the rest frame of the source, where 
9' 12 ~ 1). This reduces T 77 (e:) by a factor of r 2 ^ 1 "") where L £ « L e l ~ a at high photon 
energies (corresponding to dN p h/de oc e~ a , i.e. a is the high-energy photon index), since 
Li/ £ (i+z) needs to be replaced by Lr 2 /e(i+z) = Y 2 ^^ Li/ £ (i +z ). Second, the expression for 
the optical depth includes a factor of 1 — cos 6*12 (that represents the rate at which photons 
pass each other and have an opportunity to interact) which for a stationary source is ~ 1, 
but for a relativistic source moving toward us is ~ T~ 2 . Finally, the size of the emitting 
region can be as large as R ~ T 2 cAT/(l + z), which reduces r 77 by an additional factor of 
T~ 2 . altogether, T 77 (e) is reduced by a factor of ~ r 2 ( a+1 ), and since typically a ~ 2 — 3 this 
usually implies T > 10 2 in order to have r 77 < 1 and overcome the compactness problem. 
Using similar arguments, the lack of such a high-energy cutoff due to pair production in the 



limits on the Lorentz factor of the outflow ( 


Krolik Sz Pier 


1991 




Fenimore, Epstein & Ho 


1993; 


Woods & Loeb 


1995; 


Baring & Harding Il997; 


Lithwick & Sari 


2001). 



We note, however, that r 77 generally depends both on the radius of emission, R, and 
on the bulk Lorentz factor, T : T 77 (e) oc r~ 2a R~ l L £ a ~ 1 . Therefore, one needs to assume a 
relation between R and r in order to obtain a lower limit on the latter. Most works assume 
R ~ r 2 cAT/(l+z) fe.g- lLithwick fc Sarfc)l[ ) . which gives r 77 (e) oc r- 2 ( a+1 )(AT)- 1 L ^ 1 , 
while the lack of a high-energy cutoff up to some photon energy e implies T 77 (e) < 1. This, in 
turn, provides a lower limit on V since one can directly measure the variability time AT, the 
photon index a, and L 



4tt4(1 + z) 



^e^Fc. However, the relation R ~ T 2 cAT/(l 



does not hold for all models of the prompt GRB emission. For example, this relation does 
not hold if the prompt GRB emission is generated by relativistic magnetic reconnection 
events, with angular scales -C 1/r, that create local relativistic moti on with Lorentz factor 
7 rfi / ~ 5 — 1 relative to the average bulk value T of the emitting shell (ILyutikov fc Blandford 



20021 . |2003|). In this case AT/(1 + z) <§C R/cT 2 and the inferred value of the Lorentz factor 
from standard opacity arguments would be ~ jreiF rather than the bulk Lorentz factor of the 
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shell, T. This allows the radius of the prompt emission to be as large as R ~ 10 — 10 cm, 
close to the deceleration radius where most of the energy of the outflow is transferred to 
the swept-up external medium, and is much larger than the prompt emission radius that is 
expected in the internal shocks model, R ~ 10 13 — 10 14 cm. Therefore, we adopt a more 
model-independent approach and do not automatically make this assumption. Instead, we 
derive most of our formulas without this assumption, as well as derive expressions for T 
under this assumption, which could serve in order to test its validity. 



The Gamma-ray Large Area Space Telescope (GLAST) mission (iRita 120071 ). to be 
launched in early 2008, is expected to shed light on the high-energy emission from GRBs 
and other impulsive relativistic sources. In particular, opacity effects due to the local pho- 
ton field within the source^ are expected to be most relevant in t he GLAST L arge Area 
Telescope (LAT) energy range (20 MeV to more than 300 GeV, see iReimerl 120071 ). Thus, it 
represents a powerful tool for probing the physics of these sources. GLAST is likely to detect 
the high-energy cutoff due to pair production opacity which would actually determine T 2a R, 
rather than just provide a lower limit for it. Furthermore, in GRBs, the outflow Lorentz 



factor T may be constrained by the time of the afterglow onset (iPanaitescu &: Kumarl 12002 



Lee. Ramirez- Ruiz fc Granotll2005l ; iMolinari et al.l 120071 ). provided that the reverse shock is 
not highly relativistic, so that if GLAST detects the high-energy pair production opacity cut- 
off, the radius of emission R could be directly constrained, thus helping to test the different 
GRB models. In particular, this could directly test whether the relation R ~ T 2 cAT/ (1 + z) 
that is expected in many models indeed holds, since both R and V could be determined 
separately. This, however, requires a reliable way of identifying the observed signatures of 
opacity to pair production. This is one of the main motivations for this work. 



The leading model for the prompt emission in GRBs features internal shocks (IRees fc Meszaros 



19941 ) due to collisions between shells that are ejected from the source at ultra-relativistic 
speeds (T > 100). The shells are typically quasi-spherical, i.e. their properties do not vary a 
lot over angles < a few T" 1 around our line of sight. Under the typical physical conditions 
that are expected in the shocked shells, all electrons cool on a time scale much shorter than 



1 In the present work, we will not consider opacity effects due to interaction of high energy photons with 
the extra-galactic background light. Such an attenuation, interesting in its own right, can be added to the 
"in source" opacity in a straightforward way. Furthermore, it is expected to become significant (i.e. produce 
t 77 > 1) only at cosmological r edshifts (z > 1) and for very high photon energies (> 56 — 100 GeV at z = 1 



and > 18 — 63 GeV at z = 3; iKneiske et al.ll2004 ). and is therefore likely to significantly affect only the 



high end of the GLAST energy range, where the photon statistic might be too poor to reliably measure this 
effect. This source of opacity will be independent of time (and depends only on the redshift of the source, 
and on the photon energy), which would help in disentangling it from the time dependent opacity intrinsic 
to the source that we calculate in this work. 
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the dynamical time (i.e. the time it takes the shock to cross the shell), and most of the radia- 
tion is emitted within a very thin cooling layer just behind the shock front. Thus, our model 
which features an emitting spherical thin shell that expands outward ultra-relativistically is 
appropriate for the internal shocks model. 

As this emitting "shell" expands outward to larger radii, it builds up a photon field that 
can pair produce with high-energy photons f rom the same emission component. This effect 



has been studied in the past (see especially iBaringi 120061 . and references therein), but the 



temporal and spatial dependences of the photon field have been averaged out, corresponding 
either explicitly or implicitly to a quasi-steady state. However, in impulsive relativistic 
sources, the time scale for significant variations in the properties of the radiation field within 
the source is comparable to the total duration of the emission episode, and therefore the 
dependence of the opacity to pair production on space and time cannot be ignored, and may 
produce important effects that are suppressed in the steady-state limit. Therefore, in the 
present work we consider the full temporal and spatial dependence of the opacity, in order 
to capture all the resulting effects. 

We develop a simple, yet rich, model to investigate quantitatively the intuitive consid- 
eration that in impulsive sources it takes time to build up the (target) photon field, and thus 
the optical depth initially increases with time, so that high energy photons might be able to 
escape the source mainly at the very early part of the spike in the light curve. This results 
in a power law tail for the time-integrated spectrum at high energies, while the instanta- 
neous spectrum (which is hard to measure due to poor photon statistics) has an exponential 
cutoff. This arises since the photon energy of the exponential cutoff in the instantaneous 
spectrum decreases with time, as the opacity increases with time at all energies. Therefore, 
at sufficiently high photon energies, most of the photons escape during the short initial time 
before the optical depth increases above unity, i.e. before the cutoff energy sweeps past their 
energy. 

We perform detailed semi-analytic calculations of the optical depth to pair production, 
which improve on previous works by first calculating the photon field at each point in space 
and time, and then integrating along the trajectory of each photon. The structure of the 
paper is as follows. In § [2] we introduce our model and derive a general expression for the 
flux that reaches an observer at infinity. This expression includes the optical depth along 
the trajectory of each photon that may reach the observer, which is derived in §[3j The 
calculation of the optical depth requires the knowledge of the photon field at each point 
along the trajectory of each (test) photon. This local photon field is first expressed in terms 
of the source emissivity (§ 13.11) . Next (§ 13.21) it is conveniently rewritten as the product 
of the typical optical depth (that is approached on a dynamical time, and is similar to 



- 6- 



that derived in previous works) and dimensionless order unity expression (containing a few 
integrals) which captures the new time dependent effects that are the focus of this work. 
In § H] explicit expressions are derived for the integrands of these dimensionless order unity 
integrals. In § we derive the relevant analytic scalings for the resulting optical depths and 
observed flux, and in § [6] we present numerical results (i.e. numerically evaluate the semi- 
analytic expressions) for the opacity, light curves, and spectra (both the instantaneous and 
time-integrated spectra are addressed in §§[5] and [6]). Our conclusions are discussed in §0 



2. Calculating the Observed Flux 

2.1. Model Assumptions 

We consider an ultra-relativistic (with Lorentz factor r> 1), thin (of width <C R/T 2 
in the lab frame) spherical expanding shell, that emits over a finite range of radii, Ro < R < 
Ro + AR (i.e. the emission turns on at Ro and turns off at Ro + AR). This model can be 
associated with a single pulse or flare in the light curve. In the context of internal sh ocks 



within the outflow, AR ~ R is typically expected (IRees fc Meszarodll994j ; |Piranll2005l . and 
references therein). 

The emission is assumed to be isotropic in the co-moving frame of the emitting shell (i.e. 
the shell rest frame), and uniform over the spherical shell. In this work primed quantities 
are always measured in the co-moving frame, while unprimed quantities are evaluated either 
in the lab frame, that is the rest frame of the central source, in which the shell is spherical 
(e.g. the Lorentz factor T), or in the observer frame (e.g. the observed time and photon 
energy which suffer cosmological time dilation and redshift, respectively, relative to the lab 
frame which is at the cosmological redshift of the source). The observer is assumed to be 
located at a distance from the source that is much larger than the source size (so that the 
angle subtended by the source, as seen by the observer, is very small, and the observer can 
be considered as being at "infinity"). 

For convenience, we will use dimensionless photon energies, e, in which the observed 
photon energy, E p ^, is normalized by the electron rest energy: e = E p ^/m e c 2 . While general 
expressions will be provided when possible, we also provide detailed semi-analytical solutions 
to the model by assuming that the luminosity in the shell rest frame has a power-law depen- 
dence on rest frame photon energy e' and radius R, L' e , oc (e') 1 ~ a R b , and that the Lorentz 
factor scales as a power law with radius, T 2 oc R~ m . The approximation that V and L' £ , scale 
as power laws with radius is usually expected to hold reasonably well. For internal shocks, 
the colliding shells are expected to be in the coasting stage near the collision radius (Ro), 
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which corresponds to m = (see lPiranll2005l ; lMeszarosll2006l . and references therein). In the 
GRB afterglow, both before and after the deceleration radius, where most of the energy is 
transferred fr om the ejecta to the sho cked external medium, T (IBlandford fc McKedll976l ) 



and L' £ i (e.g., ISarilll998l : lGranotll2005l ) are expected to scale as power laws with radius. For 



GRB internal shocks, the scaling of L' e , with radius R generally depends on the details of 
the colliding shells. 

For uniform colliding shells, where the strength of the shocks going into the shells is 
constant with radius, above the peak of the uF„ spectrum, £ peak , one expects —0.5 < b < 0. 
This may be understood as follows. In this case the Lorentz factor in the shocked regions 
of the colliding shells is constant with radius, while the magnetic field scales as B' oc R~~ l . 
Therefore, since the number of emitting electrons scales linearly with radius, N e cx R, then 
L' emBX oc B'N e oc R°. The typical synchrotron photon energy scales as e' m oc B'j^ oc 
since the typical Lorentz factor of the electrons, j m , is constant for a constant shock strength. 
The energy of a photon that cools on the dynamical time (the time since the start of the 
collision) scales as e' c oc R. Therefore, above the peak of the vF v spectrum, at e' > £ pcak = 
max(e' c , e' m ), we have L' £ , = L' £ ^ mSLX (e' m /e' c )~ 1 ^ 2 (e' /e' m )' p ^ 2 oc R( 2 ~ p ^ 2 , where p is the power law 
index of the electron distribution, dN e /d^ e oc % p for 7 e > j m . Since p ~ 2 — 3 is typically 
inferred for the GRB prompt emission, this corresponds to —0.5 < b < 0. For fast cooling 



below e' 



Oc < 4) 
however, below e' k 



peak 



f' T 



(e'/e' c ) x l 2 oc R}l 2 . For slow cooling 



V > e' ) 



The simplifying assu mption of a po wer law emission spectrum [L' £ , oc (e') 1 a ], however, is 



not always valid (see, e.g.. lBaringll2006l ) . For example, in GRB internal shocks it breaks down 



for photons of energy e > r 2 /(l+^) 2 £ peak , i.e. em e c 2 > 25(l+2)- 2 (r/100) 2 (£ peak m e c 2 /100 keV)" 1 GeV. 
Indeed, photons of such energy interact with photons below the spectral break energy £break 
which is the peak of the vF v spectrum. A detailed treatment of the case of a more realistic 
spectrum for GRB internal shocks will be provided elsewhere. The exact shape of the spec- 
trum at high energies is not well constrained. Thus, we use a fiducial value of a = 2, which 
corresponds to a flat vF v (i.e. equal energy per decade in photon energy), in our detailed 
illustrative solutions, and also explore the effects of varying the value of a. 



2.2. The Equal Arrival Time Surface of Photons to the Observer (EATS-I) 



The observed normalized flux density, F e = (m P c 2 /h) F„, is calculat ed as a function of 
time and photon energy, closely following the derivation of iGranotl (120051 ). For this purpose, 
the contributions to the observed flux at any given observed time T are integrated over 
the "equal arrival time surface" (EATS-I) - the locus of points from which photons that are 
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emitted at the shell reach the observer simultaneously, at the observed time T. In the present 
work, the effects of opacity to pair production will be added at the end of this calculation, 
as detailed below. 

We consider a photon initially emitted by the shell at a lab frame time to when the 
radius of the shell is _R fj o = Rsh(to) and its Lorentz factor is I\ , at an angle of 9 tj o from our 
line of sight to the origin R = (see Fig. [1]). Due to the spherical symmetry of our model, 
there is no dependence on the azimuthal angle. The arrival time T of the photon to a distant 
observer is given by the "equal arrival time" formula: 



T 



Rt,o 



COS 9 t n 



(1) 



where the lab frame time t is related to the shell radius at that time, Rsh(t), by 

" Rsh ^ dR R sh (t) 1 f Rsb(t) dR 



(3c 



2c 



T 2 (R) 



+ o(r- 4 ) . 



(2) 



In Eq. ([T]), T = is chosen to correspond to a photon that is emitted at the origin at to = 0. 
Eq. ([2]) relates t and R s h(t), so that the locus of points (R t fi,@t,o) that keep T constant 
defines the EATS-I at time T. For a coasting shell (m = 0), it is a well-known result that 
the EATS-I is an ellipse^] of semi-major to semi-minor axis ratio T (Regs 1966 ). The flux 
density at the rescaled energy e is obtained by integrating over the luminosity in the shell 
rest frame, L' £ ,, along the EATS-I JGranotlbopj ): 



(1+*) 

Andl 



5 3 dL' e , 



(1 + *) 
8vrd| 



dy^5 3 (y)L' £/ (y) 
ay 



(3) 



where 8 = (1 + z)e/e' is the Doppler factor of the emitted photon (between the co-moving 
and lab frames), /i i)0 = cos 6^0 is the cosine of its angle of emission, and we defined the 
normalized radius y = R tt o/RL, where Rl = Rl(T) is the largest radius on the EATS-I at 
time T. The integration is performed along the EATS-I, and the boundaries for y are 



Vmin{T) = min 



Ro 



Rl(T)\ 



mm 



R + AR 
Rl(T) _ 



(4) 



since the emission turns on at Rq and turns off at -Ro + AR. For the times T relevant to 
the problem, corresponding to the arrival of photons to the observer, Rq/R l (T) is always 
smaller than 1. 



2 It actually represents an ellipsoid, keeping in mind the symmetry around the line of sight to the center 
of the emitting spherical shell, and the lack of dependence on the azimuthal angle. 
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In order to evaluate the integral above, we now derive expressions for the integrand. 
Defining r L = T(R L ), T 2 oc R~ m can be rewritten as T 2 (R)R m = T 2 L R™ = constant, and 
thus T 2 = T 2 L y~ m . Eq. (|2J) now becomes 



, _ Rt,o Rl y m+1 n( V -^ 



(5) 



In the limit of small angles (6^ <C 1, which is relevant for T ^> 1), Eq. ([I]) implies tn—Rt,o/c 
T/(l + z) - Rt fl 9 2 /2c, which together with Eq. © yields 



T 



RlV 



m+l 



+ 



RtflQlo 



(I + 2) 2(m+l)r|c 2c 

As can be seen in Fig. [TJ a photon that is emitted at R t $ = Rl (corresponding to y 
Rt,ol Rl(T) = 1) remains along the line of sight (9 t = 9 t) o = 0), so that Eq. (jSJ) yields 



(6) 



R L (T) = 2(m+l)T 2 L [T/(l + z) 



cT 



l + z) 



where T = T(Rq), and can be rewritten as 



Ro 



-1 



Tn 



l/(m+l) 



Tn 



(1 + z)R n 
2(m + l)cT 



> • (7) 



(m + 1)T 2 



We have introduced the time T at which the first photons reach the observer (corresponding 
to a photon emitted at Ro along the line of sight, 9 = 0): Rl(T ) = Rq. Since fj, tt o ~ 1— o/2, 
Eq. (jSJ) implies 



cfy 2(m + l)r| 
Finally, the Doppler factor of the emitted electron is given by 

1 _ 2r _ 2(m + l)Y L y- m ' 2 

= r(i - (3 cos ^,0) w 1 + (r^ ) 2 ~ m + 1 ' 

and its value at .R^ (which corresponds to y — 1) is 5(Rl) = 2T L . Since 



(9) 



(10) 



" e' " 


l-a 











'IF 



where e' = (1 + z)e/5, we obtain: 
L £ > = L( 1+z ^ E / 2 r L (RL 



5 



2T, 



y b = l\ 



(l+z)e/2r (-^o 



2r ; 



a-l 



Rl 
Ro 



b-m(a- 1)/2 



(12) 
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The effect of pair production opacity will be treated in this work in a somewhat simpli- 
fied manner, by assuming that photons which pair produce do not reach the observer, and 
ignoring the additional opacity that is produced by the secondary pairs and the photons 
emitted by these pairs. Under these simplifications, the effects of opacity to pair production 
can be included by adding a term exp(— r 77 ) into the integrand in Eq. ([3]), where T^yy IS 3j 
function of y, e, AR/R , and T/T . Thus, by combining eqs. ([9]— [12]) with Eq. <^fy, we obtain: 

{1 + z) fT Y b - maWim+1)] 

x f ^ dy ( m + \ \ y b -^ e-^ , (13) 

where Eq. (J7J) is used to derive the scaling Rl(T)/R q = (T/T ) 1//(rra+1 ), and 

/ AR T L \ 

r 77 — r 77 I y, c, , , 2a I , 
V -^o 1 o 1 o -"-o / 

as is shown later on, where To = T(Ro), and L £ « Lo£ 1_a is the observed isotropic equivalent 
luminosity. Unless specified otherwise, the derivations throughout this work are valid for a 
general value of m. For a coasting shell (m = 0), which is a case of special interest (as it is 
expected, e.g., for internal shocks), Eq. ([TBI simplifies to 



(1 + 2) ( T \ b /" ymax 
F e (T) = (^) / «*» 

\ w / u ymin 



a+b e -r 77 



We have expressed the observed flux density for our model as a function of the observed 
time T, and we now need to derive the expression of the optical depth r 77 . We gather here 
the dependence on y of two quantities that will be needed later on: 

Ro = -B~ = = —— = - — , X = (r t , O tjO ) = , — • 15 

i2t,o y AR Rt,o y\T J (m+1) 

In order to facilitate reading, we include in Table[TJthe most common quantities used through- 
out this work. 



3. Computation of the optical depth 



As in the previous section, we consider a "test" photon emitted by the shell at radius 
R t fl and angle 9 t<Q with respect to the line of sight (see Fig. [TJ). All the quantities with a 
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subscript T will always refer to such a test photon. We wish to calculate its optical depth to 
pair production with all the other photons which are emitted by the same source and denoted 
by a subscript 'i' (for p otentially "int eracting"). The differential of the optical depth to pair 



production is given by (I Weaver! 1 19761 ) 



dvt ' 

dr 77 = a*[x(e t ,e h Li t i)](l - n ti ) ^ ^ dVLidEids . (16) 

In this equation, ds is the differential of the path length along the trajectory of the test 
photon; n,- n and Ei = £im e c 2 are the number density, solid angle, and photon energy 
of the photon field along the path of the test photon with which it might interact. B For 
convenience, St and Ei denote the values of the corresponding dimensionless photon energies 
in the lab frame, rather than in the observer frame (as is the case for e), i.e. without the 
cosmological redshift, so that s t = (1 + z)e should eventually be used in order to evaluate 
the optical depth at an observed value of e. The Lorentz invariant cross section for pair 
production cr*(x) is 

°*(X) = fe 4 + 2x - 1)HX + Vx^) - X(l + X 2 ) Vx^l] , (17) 
X b L J 

X = "" ) ■ (18) 

where \ is the center of momentum energy (in units of m e c 2 , of each particle - each of the 
two interacting photons, and the produced electron and positron), and \i t i = h t ■ hi is the 
cosine of the angle between the directions of motion of the test photon (h t ) and a potentially 
interacting photon (hi). In order to evaluate /i ti , we need to specify the geometry for our 
model: a spherical emitting shell, whose emission depends only on its radius _R s h (i.e. at 
any given radius its local emission does not depend on the location within the shell) and 
is isotropic in its own rest-frame. Under these assumptions, the radiation field will depend 
only on the radius R and the (lab frame) time t, and at any given place and time it will 
be symmetric around the radial direction (see Fig. [2]). Therefore, at any point along the 
trajectory of the test photon, we can use a local coordinate system, S r , whose z-axis is 
aligned with the radial direction (from the center of the shell to that point), z r , and such 
that the direction of motion of the test photon is in the x-z plane. In this frame the polar 



3 We do not add a facto r of 1/2 due to double counting (as was done by, e.g., iBaring fc Hardinall997 ; 
Dernier fc Schlickeiserlll994l) , as it should not appear in the expression for the optical depth. We discuss this 
point in more details in annex lEl 



angles are denoted by (9 r ,<f) r ), and 

h t = x r sin 9 t + z r cos 9 t , (19) 
hi = x r sin 9 r cos 4> r + y r sin 9 r sin r + z r cos 6V , (20) 
[Hi — hi • h t — sin 6*t sin 9 r cos r + cos 9 t cos 6* r . (21) 

Note that 9 t varies only with s. The integration over the solid angle in the lab frame in 
Eq. (1161) can conveniently use the frame S r which is at rest in the lab frame, i.e. dQi = 
dQ r = d<f) r dfi r . The optical depth of the test photon is then given by: 

T^j(e t ,9 t ,o,Rt,o) = J ' ds J ' dei J dQ r a*[x(e t , e h /jl u )] (1 - Hu) d ^^ £ , ■ ( 22 ) 

Next, we express the derivative in the integrand of Eq. fl22|) . which represents the photon 
field along the trajectory of the test photon, in terms of the source emissivity. In addition, 
we make a series of changes of variable in order to simplify the expression for the optical 
depth. 



3.1. Expressing the photon field in terms of the source emissivity 

In § I2.2[ we expressed the observed flux as an integral over the EATS-I of photons to the 
observer at an observed time T. These photons travel along straight line trajectories that pass 
through the photon field. As a result, we integrate the contribution to the optical depth at 
each point along the path of each photon, treating it as a test photon. This is the integration 
over ds in Eq. (122]) which, as we show below, can be replaced by an integration over dR t . 
In the other two inner integrations Rt is kept fixed, and the photon field, drii/dQ r dei, needs 
to be evaluated as a function of £j, [i r and R t . For a given test photon that is emitted at 
(i?f , Ht,o), the value of R t also determines the value of the lab frame time t t . We remind 
the reader that Rt and t t are always computed in the lab frame, and that Rt is in general 
different than R s h(t t ), i.e. at a general time the position of the test photon does not coincide 
with that of the shell. We proceed first to relate the photon field at (t t , Rt) to the emissivity 
in the local frame of the emitting shell, which is easier to specify, and simpler. The Doppler 
factor of the emitted photon is given by 

^7 = rn * S =r(i + /3/4), 

e\ T(l - Pud 

where /ij = cos 9i — /3 • hi and = cos 9[ = (3 ■ h\ are the cosines of the angle between the 
bulk velocity of the emitting fluid (f3) and the direction of the interacting photon in the lab 



(23) 



-13- 



frame (hi) and in the comoving frame of the emitting fluid (nQ, respectively. Furthermore, 

, = Hi-fi ^ dty = dty = dfj/_ = §2 ^ 
1 1 — P/ii dVl r dVti dji 

since dQi = dfyd/Ai and = fa. We are interested in the differential density of photons of 
energy and direction of motion in the solid angle dQ r around the direction hi, which is at 
an angle 9 r from the radial direction, at a radius R t and time t t . This density is related to 
the specific intensity of the photon field by: 

^ " dsdtlm, = £ * mec3 dBk {fli) ' (25) 

where the (normalized) specific intensity I £i is the energy (dE) per unit normal area (dS 
where dS/dS = h), per unit time (dt), per unit (normalized) photon energy (dei), per solid 
angle (dfli = dfl r ) around some direction hi of the (potentially interacting) photons. 

The differential (normalized) specific luminosity (in our case, from a small part of the 
emitting shell) is defined as dL £ = dE/deidt, while the isotropic equivalent (normalized) 
specific luminosity is defined by: 

dL £iMo = 4tt^ . (26) 

The contribution of an emitting element with dL £uiso to the (normalized) flux density dF £ . = 
dE / dSdtdei and to the (normalized) specific intensity I £i at a point located at a distance r 
from it is 

dF £i = = I £i (h)dQ r , (27) 

and is along the direction h from the emitting element to that point (i.e. here dS is the 
differential of the area normal to n, dS = h-dS). Finally, we can conveniently express dL £ ^ iso 
in the comoving frame (i.e. the local rest frame of the emitting shell), 

dL F . dE , dE 1 , ^ , ,„„ N 

= te ii = 4 "ffl- = * 4l ;™; = SdL " • (28) 

where the last equality follows from the assumption that the emission is isotropic in the 
comoving frame. Because the emission is assumed to be uniform throughout the shell, dL' e , 
depends only on the radius of emission of the potentially interacting photon, R e , and not on 
the location within the shell. Apart from the emission radius, R e , the position of an emitting 
point on the shell is also specified by the polar angle, 9 e , which for convenience is measured 
with respect to the direction from the center of the sphere to the location of the test photon 
(at a radius R t ) where the flux (or some other property of the photon field) is calculated 
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(see Fig. [2]). As a result, we can write dL' £ , = L' s ,(R e )dp, e /2, where \i e = cos9 e and L' s ,(R e ) 
has been defined and discussed in §[2l We finally combine eqs. (j25il . f |27l) and ( |28l) to obtain 
the expressions for the normalized specific intensity, 



47T Atit 2 



djj, r 



(29) 



and the expression for the photon field which appears in the integrand in Eq. fl22|) . 



drii 
dejdQ r 



(4n) 2 em e c 3 r 2 



dfi e 



dfi T 



(30) 



The derivative in the last term to the right of these equations must be computed along the 
equal arrival time surface (EATS-II) of photons to R t at t t , where I £ or driij 'deidVL r are to 
be calculated. We can now rewrite Eq. fl22|) as: 



9t,o, Rt,o) 



(47r) 2 m e c 3 



jdsjdejdn ^ £ ^^\ l-, tl 



dfX e 



dfi r 



(31) 



We have thus replaced the photon field by the specific emissivity in the expression for the 
optical depth. The boundaries of integration will be specified explicitly later on. We now 
want to simplify this triple integration in order to make it easier to evaluate. 



3.2. Analytical reduction 

In the remainder of this work, we will make use of various dimensionless radii, which 
are gathered in Table [TJ and greatly simplify the analysis. Furthermore, it is much more 
convenient to work with such quantities of order unity inside the integrand. We thus rescale 
R e and R t by introducing R e = R e /Rt and R t = Rt/Rt,o- Furthermore, the notations 
R = R/R t and R = R/Rt,o will be used for other rescaled dimensionless radii as well. While 
clearly 1 < R t < oo, the range of R e is much more complex and will be extensively discussed 
in §|H For now, we want to simplify Eq. (13"TT) by changing integration variables. We give 
here the main results and leave the details of the derivations for Annex [A] 

As has been mentioned above, the integration over ds can be replaced by an integration 
over R t . Under the approximation of large Lorentz factors (r ^> 1), and thus small emission 
angles (0 tiO 1), one obtains ds = Rt^dRt (see the discussion following Eq. [A2] for more 
details). Besides, since we integrate over dQ r = d<p r dfi r and the integrand contains \d[i e jd[i r \ 
we can conveniently change the integration over fi r to an integration over R e . We show in 
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the Annex lAl that 



dfl e 



dfX r 



dfj, r = dR e 
dRp 



since dfi e /dR e > 0, where the limit of integration over R e should be in increasing order (i.e. 
the integration should be from small to large values of R e ). The optical depth now reads: 



ax 



00 dR t r °° 



(An) 2 m e c 3 R tfi J 1 R 2 j 2/e 

x 



dei 



2tt 



dR P 



Ro/Rt 



Hi) 



Si f 2 dR 



(32) 



Next, we can follow the hind-sights of IStepney fc Guilbertl (119831 ) and iBarind (119941 ) in 
order to cast the integrations over (d<p r , dei) into a much more practical form. In order to 
perform this change of variables, it is necessary to specialize the specific luminosity to the 
dependence discussed in §H L'(R e ) = L' Q h(R e l RoMei') 1 ^ = T Q a L Q (ei') l ~ a x h(R c R t /R ), 

i 

where h is a general function of R £ /Rq that satisfies h(l) = 1 (for details see Appendix IA.3I) 
and L' = L' £ , =1 (Rq). Note that L = TqL' is approximately the observed isotropic equivalent 
luminosity at a photon energy of m e c 2 ^ 511 keV, near the peak of the spike in the light 
curve which corresponds to the emission episode that we model for AR ~ Rq. 

For convenience, we rescale all the quantities in the integrand of Eq. fT32|) which are 
not of order unity by the relevant power of the Lorentz factor at radius R t , T t = T(R t ), so 
that the rescaled quantities (which are denoted by a bar) will be of order unity. We rescale 
5 = S/T t and djx e = T 2 d[i e , but do not rescale f which is already of order unity. Thus, 



\2+a 



d[i e 
dR P 




ma/2 — 

6 2+a 



djji e 
dR P 



(33) 



and the expression for the optical depth becomes: 



r 77 (s t , e t , , Rw) = r,er l Rl~ ma ' 2 I [ dRj-^- ^ h ( R*^ ] C^(C) , (34) 



1 R, 



° f 2 dR e 



where 



(r5°Bo,a,L„) 



T* 0.402 (|r V3 10^") L " 



487r 3 m P c 3 a 5 / 3 K 







',2 



(To,2) 2a R 



0,13 



(35) 



and Lo, 52 = L /(10 52 ergs" 1 ), R , 13 = i? /(10 13 cm), r 0)2 = r /100, C = (C+ ~ C-VC-, 
( + = [1 — cos {0 r + t )} /2, (- = [1 — cos (9 r — 9 t )] /2, and H a (z) is a function discussed in 
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Annex [AT3l In Eq. fl35l) . r*, the only quantity requiring astrophysical input, is a constant of 
the order of the optical depth to pair production at a photon energy of m e c 2 at R in quasi- 
steady state (near the peak of the spike in the light curve for AR ~ R ). Note that since 
both the photon index a and L (roughly the isotropic equivalent luminosity) are observable 
quantities (the latter requiring knowledge of the source redshift), the observation of a high- 
energy spectral cutoff due to pair production opacity can enable the determination of Tl a R . 
In the limit of small angles that is appropriate for large Lorentz factors, is of order T~ 2 , 
so we define (L = rf£_, where T t = T(R t ) = T R^ m ^ 2 . Thus, 

(C-)° • (36) 




\b 



Under the assumption that h is also a power-law of index b, h(R e ) = (i? e /i? ) 
(R e R t /Ro) b , the expression for the optical depth in our model simplifies to: 

r 77 (£ t , 9 t ,o, Rt,o) = T (e t , Rt,o)T{x) , (37) 
T (e t ,R tfi ) = T*er 1 R 1 - b - ma/2 , (38) 

Hx) = dR t ^" 2 +"W 2 [ d R e 5 -^l . fajj* CH a (C) . (39) 
J i J r dR e 

In order to proceed further, we need to obtain an explicit expression for the innermost 
integrand of J 7 , by a detailed examination of the geometry of the photon field. The next 
section will be devoted to this analysis, which constitutes the main novelty of this work. We 
will evaluate the optical depth (Eq. [2Z])? taking into account that the photon field is not 
homogeneous along the test photon trajectory, but the contribution to the photon field is 
actually built up in time. 



4. Calculating the Photon Field 

4.1. Equal Arrival Time Surface of Photons to the Test Photon (EATS-II) 

In this section we calculate the photon field at a general radius Rt and time t t , along the 
trajectory of a test photon. For this purpose we need to consider the contribution from all 
photons that arrive at the instantaneous location of the test photon, (R t , t t ), simultaneously. 
The locus of points where all such photons are emitted, taking into account that the emission 
occurs only in the shell, forms a two dimensional surface referred to as the equal arrival time 
surface (EATS-II) of photons to the instantaneous location of the test photon. The local 
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photon field at (R t , t t ) is calculated by integrating the contributions over this surface. We 
stress that this surface (EATS-II), is different from the equal arrival time surface of photons 
to the observer at infinity (EATS-I). 

Fig. [1] shows the basic configuration for our calculations and illustrates the relation 
between the two different equal arrival time surfaces (EATS) of photons: 1. to the observer at 
infinity (EATS-I), 2. to the instantaneous location of a test photon (EATS-II). It can be seen 
that the EATS-II grows with the lab frame time t, and therefore also with the radius of the 
test photon R t . Furthermore, each EATS-II encompasses all other EATS-II corresponding 
to smaller times, and is encompassed within all the EATS-II which correspond to larger 
times. In particular, all EATS-II are within the EATS-I, which corresponds to the limit of 
the EATS-II for an infinite time (when the test photon reaches the observer at infinity). All 
of the EATS-II and EATS-I pass through the emission point of the test photon, and for case 
2 and 3, also through the place where the photon crosses the shell (i.e. its location in case 
2). These are general properties of the EATS-II. 

We now proceed to calculate the EATS-II and the expressions for relevant quantities 
along this surface, which are needed in order to calculate the local radiation field. From the 
geometry of our problem (see Fig. [2]), we can immediately derive the two following equations: 



where R e = R e /Rt and f = r/R t . The equal arrival time surface (EATS-II) of photons to 
(R t , t t ) is determined by the condition that r = c(t t — t e ) = c[t t — t s h(-R e )], where the photons 
are emitted at a previous time t e when the shell is at a radius R e = i? s h(t e ). The EATS-II 
equation is thus given by 



which relates the radius (R e ) and angle (8 e = arccos fi e ) of emission along this surface. 

The expression for t sh (R e ) depends on our assumption about the expansion of the shell. 
If the latter occurs at constant speed, then t s h(-R e ) = R e //3c and in the limit of R t — > oo 
Eq. ( j4"2l reduces to j3(ct t — Rt) = -R e (l — flUe), which is the usual polar equation of an ellipse 
(setting cT = ct t — Rt). In this simple case, using the short notation R sh = R s h(t t ), we have 
f = {Rsh — R e )/P, and the EATS-II is given by 



r 2 = 1 + Rl - 2R e n e = (1 - R e f + 2R e {\ - Me ) , 
R\ = 1 + f 2 - 2rn r , 



(40) 
(41) 




(42) 




(43) 
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while the lower and upper limits for the range of R e values along the EATS-II, which corre- 
spond to n e — — 1 and \i e = 1, respectively, are given by 

(R sh + (3)/(l + (3) R sh >l 



1 + 



(44) 



,(it! S h-/3)/(l-/3) i?sh < 1 • 

Note that we have not assumed T = (1 — /3 2 ) -1 / 2 ^> 1, so these results are valid for an 
arbitrary velocity, as long as it is constant with radius. 

Combining eqs. (T40]) and (jUJ) we also obtain 



//,. 



-R e /i e 



c 



1 — R e fi e 



(45) 



[l-R e y + 2R e (l-fi e ) 

where in the last equality we have also used Eq. (I42p . so that it is valid only along the 
EATS-II (while the first equality is valid more generally, as it is derived directly from the 
geometrical setup). 



4.2. Radial Dependence of Relevant Angles, /i e (-R e ) and /i r (_R e ), along EATS-II 



Specifying for 1 -C T 2 = T 2 i? m , we can rewrite Eq. ([5]) as 

R m+1 



t S h(R) = — 

c 



R + 



2(ra + l)Tf 



+ o (it 4 ) . 



(46) 



Thus, Eq. ( 1421) implies 

f = R sh - R e + 



pm+1 Dm+1 

2(m + l)r 2 



+ o (it 4 ) 



f 2 



R s h — Re) + 



+ o (it 4 ) 



(47) 



(4* 



(m + i)r 2 

Note that i? e < i? s h, because i e < t 4 (due to causality) and R s h(t) is an increasing function 
of t. The equality only holds when R t = i? s h(tt), i.e. when i? s h = 1 (case 2 below). Thus, 
eqs. (1401) and (1481) give (to the order of IT 2 ), 



2T 2 t (l-fi e ) = (T t e e ) 2 
1 

R P 



i? s h — -R, 



1 -R f 



Rsh — Re ) ( -R 



L sh 



(m + 1) 



(49) 
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The two terms on the right hand side of the equation are typically of the same order since 
| i? s h — 1| < a few T^ 2 , i.e. R s h = 1 R s h(U) = Rt- This immediately implies 



1 



dR e ~ 2Y\Rl 



R 



sh 



m + 1 



RT 1 + 



m+l 



- R 



m+2 



(50) 



Now we turn to \i r . From eqs. (jHJ) and (flBl) we obtain 



-Rsh — R e I + 1 1 — R\ 



2 R 



sh 



+ 



pm+l pm+1 

4(m + l)rf 



1 — 



o (17 4 ) 



(51) 



dR P 



dr 



dR P 



l-R 2 - r 2 



2r 



Re 

r 



r: 



2[R sh -R e ^ 4r * 



l-Rl 



R s h — R t 



- 1 



+ 



pm+1 pm+1 \ / p p I 



2(m + l)r t 2 [R sh -R ( 



+ ^ (rr 4 ) 



where 



dr 



pm 



1 + 



1 



1 



(52) 



(53) 



di? e ' 2r 2 ' 2T 2 (R e ) (3{R e , 

This can easily be understood since r = c(t t — t e ) along the equal arrival time surface, so 
that dr = —cdt e and df/dR e = dr/dR e = —cdt e /dR e = —cj (dR e /dt e ) = —l/(3(R e ). 

The maximal radius of emission, R e ,ma,xi from which a photon reaches a point at radius 
R t at the time t t is determined by the photon that is emitted at 9 e = (i.e. /i e = 1), along 
the line connecting that point to the center of the sphere. Thus, 



1 -FL 



-Rsh — Re 



C 

Rt 



tt — tsh(Re,r 



+ 



pm+1 pm+1 



2(m + l)r 2 

and the problem naturally divides into three cases. 



+ o (r- 4 ) 



(54) 
(55) 
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4.3. Properties of EATS-II According to Relative Location of Test Photon 

and Shell 

The properties of the EATS-II qualitatively change according to the location of the test 
photon relative to the shell at the same lab frame time, t t . Thus the problem naturally 
divides into three illustrated in Figs. Q] and El If the photon is emitted at an anglej 

9tfl > l/Ttp, i.e. x = (TtflOtfi) 2 > 1, it initially lags behind the shell (case 1), since due 
to the aberration of light (also referred to as relativistic beaming) this corresponds to an 
angle greater than 90° from the radial direction in the co-moving frame of the shell. The 
photon eventually catches-up with the shell and crosses it (case 2), since the latter is moving 
at a velocity slight smaller than the speed of light. After it crosses the shell, it remains 
ahead of the shell (case 3). A photon that is emitted at 6 t ,o < l/r t) o, corresponding to 
x < 1, immediately gets ahead of the shell (case 3). All photons are always emitted at the 
shell, so the point of emission is considered case 2. Like the later shell crossing for photons 
with x > 1, case 2 corresponds to a single point along that trajectory of the test photon, 
unlike cases 1 corresponds to a finite path along the trajectory, and case 3 corresponds to a 
(practically) semi-infinite interval (as far as the observer is considered to be at "infinity" ; the 
contribution to the opacity at large distances from the source, however, becomes negligible). 
The three different cases are discussed in detail below, and the relevant expressions for each 
case are derived. We start by defining some useful quantities for this purpose, which will be 
very helpful later on. 

In the limit of small angles, Eq. flAll) yields 

{YAf « xR- t m - 2 , (56) 

where x = {T t fi6t.o) 2 is the square of the normalized emission angle of the test photon. 
Evaluating Eq. (JIB) at = R s h(t t ) gives 

t = ^ + nSwi + ° (r * 4) • (57) 

which can be rewritten in terms of the quantity 

f m = 2(m + l)T 2 t 0| - l\ = 2(m + l)V 2 t (R sh - 1) + + O (Y^) , (58) 
that plays a major role in the following derivations. 



4 More generally, the condition is cos Otfi < (3, but for T t ,o ^S> 1 and 9 t ,o "C 1 this reduces to # tj o > l/T{,o- 
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For an emission episode starting at Rq = 0, the inequality ct t > R t is required in order 
to have a non- vanishing radiation field at the point (R t ,t t ). If the emission turns on at a 
non-zero radius Rq, this condition generalizes to 



C^i _ j y Rq _ R™ +1 

Rt ~ 2(m+l)r 2 (i? ) ~ 2(m+l)r? ' 



17 ~ 1 > o/ , v = 0/ " tt^ • (59) 



This implies that f m >0 (for m > — 1, which is assumed in this work, and is typically the 
case for the astrophysical sources of interest). 

We note that f m < 1 for R^ < 1 (when the test photon is traveling in front of the shell), 
f m > 1 for R sh > 1 (when the test photon is traveling behind the shell), and f m — 1 for 
i? s h = 1 (when the test photon is at the shell). It is convenient to express f m as a function 
of our primary variables. Using 

R\ = Rl + z 2 = R 2 tfi sin 2 6 tfi + [Rtjo cos 6 i)0 + c(t t - t )} 2 , (60) 

where R± is the distance between the line of sight to the origin and the trajectory of the 
test photon (see Fig. [I] and Eq. [Al] ) and solving this second order equation, one obtains 



c(t t -t ) (h \ / , ft 



v- 

7 f,0 



R t,o v 1 \ 2Rt. 
Recalling that ct /R tfi = 1 + l/2(m + l)r 2 + 0(T^), we finally obtain: 

/ rf \ 1 + x(m + 1) f 1 - A -1 
/ m (it* t ) ^ 2(m + l)r 2 - 1 J = L + (r- 4 ) . (62) 

Fig. H] shows the dependence of f m (Rt) on the parameter x = (r tj0 #t,o) 2 - For fi f = 1 we 
always have f m = 1 since the test photon is emitted at the shell. For x > 1 the photon 
initially lags behind the shell (case 1), and the equation f m — 1 that can be expressed as 
R™ +2 — [1 + (m + l)sc] i?t + (m + l)a; = has an additional non-trivial solution, R 2 , which 
corresponds the the point where the photon crosses the shell. For m = and m — 1, it is 
given by R 2 = x and R 2 = (y/1 + 8x — l)/2, respectively. 



4-3.1. Case 1: Test Photon Behind the Shell, R t < R s h{tt) 

In this case 

, . , / a few\ ~ ~ . , „ a few . , 

Rt < Re^ex < Rsh(tt) < Rt ( 1 + -p2-J ^ 1 < #e,max < #sh(*t) < 1 + -Jq~ , (63) 
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where the last approximate inequality holds for emission angles (r t 6 e ) 2 < a few, from which 
most of the contribution to the observed flux arises, and are therefore the ones of relevance. 
An expression for R s h{tt) may readily be obtained through (see Eq. [16] ) 



ct 



R s h_(t 



ct, 1 



*-<«*) = £ 1 - + ° ( r r 4 ) = £ - ,7-^ + o (r-) , (64) 



Rt 2(m+l)rf v * ' R t 2(m+l)r 

while -R ei max is obtained by equating the two expressions for f, from Eq. (HTj) and Eq. (HOl) 
for /i e = 1, 

pm+l pm+1 



mm - 1 l e,max 



-Re, max 1 — -Rsh -Re, max "I - „/ , _, n-,-,9 "I - (-TV ) — -Rsh -Re, max "I - (-TV ) 

2(m + 1)1 1 



(65) 



which implies 



2(4, m »-l)«(^-l)«g-l)- i5 ^ ? , (66) 

- = Kt + - iraff + ° ^ - H 1 ■ w 

While 6 e is always small, (T t 6 e ) 2 < a few, in the case studied in this subsection 6 r can 
range from zero to 7r and it is not obvious a priori whether it can be taken to be either 
large, (r^,,) 2 ^> 1, or small, (T t 6 r ) 2 < a few. We argue that when 6 r is large, the photons 
must be emitted at a large angle relative to the direction of motion of the emitting shell 
{6i = 9 r + 9 e ), and are therefore significantly suppressed by relativistic beaming. This effect 
wins over the increase in the reaction rate due to the larger angle between the test photon 
and the interacting photons, that is manifested by the factor of (1 — fi t i) in the integrand 
for the optical depth. Therefore, the dominant contribution to the optical depth occur from 
small 8 r values, and we can therefore make the approximations that are appropriate for 
(r t r ) ^ a few. We express these considerations more quantitatively in Annex [B] Thus, we 
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obtain: 



(r^ e ) 2 = 2r t 2 (i-^, 
'i-Re 



-Re, max R( 



R 



m+l pm+1 



(m + l)R e 



(m + l)R e 

f m (Rt) - R? +1 } + o (r- 2 ) 



R™ +1 +4(m+l)T 2 t i? e , max -l 



(r^ r ) 2 = 2r 2 (i- /Ur 



- + 2(m + l)r 2 ( i? sh - 1 



iL 



fm(Rt) — R, 



m+l 

e 



(m + 1) ( 1 - R 



(to + 1) ( i? sh - iL 

o (r- 2 ) , 



+ ^ (rr 2 ) 



^ (rr 2 ) 

(68) 



(69) 



d/ie _ 1 

dR e ~ 2(m + i)r 2 I ^ 



m+l 
e,max 



iC +1 +4(m + l)r 2 (iL, max -l 



+ (m + l) J R™- 1 LR, 



R f 



f m (R t )-R™ +1 + (m+l)R 



r +1 ( 1 - R f 



2(m + l)r t 2 i?2 



d (m + l)iL™ +1 ( i? sh - R e j - R sh R2 +1 - RT +1 + 2(m + l)T 2 t (R sh - 1 



dR P 



2(m + l)r 2 (# sh - R t 

[m + 1)^ +1 ( 1 - Re ) - fm(Rt) + R? +1 



2(m+l)T 2 t (l-R 



(70) 



(71) 



max ) = 1, since ( 1 - R e 



We note that, as expected, jJ. e {R e , 
djji e /dR e > 0. The Doppler factor is given by 



(l - i? e ) = (\ max - R e ) +0 (IT 2 ) while 



2T 



2(m + l)T t R 



(m+2)/2 



l-R, 



1 + T 2 (9 e + 6 r ) 2 (m + ^rn+l (l _ + f m (fy _ £m+l 



(72) 



where we have used eqs. (168]) and (1691 as well as T 2 = TfR e m ^> 1 and # e + r <C 1. Finally, 
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r « 1 — i? e , and thus 

^ a+2 ^ e _ g"+ 2 ^ 2(2r f ) a (m+ l) 1+ ^e a+f(2+Q) (l -R c ) a 

rl ' dR e * fl ' dR e ~ r ~ ' ~ ^ - ~ ^ 1+Q ' l " Jj 



(m + (l - 4) + /m(4) - ^ e m+1 



4-3.2. Case 2: Test photon at the shell, R t = R s h{tt) 

This is a limiting case between case 1 and case 3, when the test photon is located on 
the shell: t t = t Bh (R e;max ), f min = 0, and # e , max = 1, i.e. 

Re,max = Rsh(U) = Rt , -R e ,max = Rsh(tt) = 1 • (74) 

This means that the last emitted photons that still reach the point (Rt,tt) are emitted at 
that same point in space and time, i.e. the equal arrival time surface ends at that point. 
Therefore, 



1 - R e 1 - R™ +1 

{Y t 9 e f = 211(1 - fi e ) = ^ '- , (75) 

(m + l)R e 



R e [l-R 



■m+l 



(TAY = 211(1 - f^) = ^— ^ , (76) 

(m + 1) ( 1 - R e 



dfi e 



mRT +1 1 - R P + 1 - RT +2 



dR e 2{m + l)T 2 t R 2 e 



dR P 



2(m + l)T 2 t (l-R e y 
In the limit where iL ~ 1 (i.e. 1 — iL <C 1) we have: 



(77) 



(7* 



,2 /-, r) \ 2 /-n /i \2 _ 1 m + 2 / ~ \ d/i e 1 — i? e 9 e dfi r (m + 2) 



(T t e e y mil -Re) , (r^r « 1 — 1 - R e . . . - „ ~ . - ~ , . 

(79) 

In this limit f ~ 1 — i? e , which implies [see Eq. (1301) ] that J e oc \dfi e /dfi r \/f 2 oc (1 — -R e ) _1 , 
i.e. the specific intensity diverges at the angle 6 r = 6 r>raax = 1/Y t , and vanishes above this 
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angle. This can be understood as follows. In this limit f < 1, i.e. r <C R t = R s h{tt) and the 
curvature of the shock front becomes unimportant, so that in order for a photon to reach the 
point (R t ,t t ) together with the shock front it must propagate along the shock front, which 
corresponds locally to an angle of 1/T (or more generally cos 9 = (3) from the normal to the 
shock front, i.e. the radial direction in our case. 



4.3.3. Case 3: Test Photon Ahead of the Shell, R t > R s h(U) 
With the causality condition -R e ,max < Rshitt) and Eq. ( |59l) . we now have: 

Rq < -Re,max < Rsh(t t ) < Rt < Ct t 



RtK +1 



2(m + l)r t 2 ' 



>R < -R e ,max < Rsh < 1 < 



Ctt 



pm+1 



R t 2(m + l)r? 



As a result, Eq. (Jool) yields: 



1 -^?e, max Rsh Re 



+ 



pm+1 pm+1 | 

sh - n 'e,max I 

2(m + 1)T? 



+ o (r- 4 ) 



1 - R. 



sh 



pm+1 _ pm+1 pm+1 _ pm+1 

"sh -"^max ^ -"-sh - rt '0 ^ 



pm+1 



2(m + l)r 2 



< 



1 



2(m + l)r 2 " 2(m + l)Tf 2(m + 1)T 



2 ■ 



and R P 



R 



sh 



2(m + l)r, 2 



2(m+ l)r 



pm+1 



2 / °tt 
1 1 



1 — -Rsh 

l/(m+l) 



l/(m+l) 



1 



fm(Rt 



l/(m+l) 



(80) 
(81) 



(82) 
(83) 



(84) 



Taking these results into account, we now derive the relevant expressions from eqs. ()49l— 1521) . 
The leading terms for 1 — /x e , 1 — /i r , and their derivatives with respect to R e are all of the 
order of O (r^ 2 ) . Thus, we obtain: 
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(ry? e ) 2 = 2r 2 (i-^ e ) 



(ry? r ) 2 = 2r 2 (i - /i r 



(l - RA (i?™+ix - RT +1 ^j 



1-R f 



(m + l)R e 

fm(Rt) — R™ +1 



+ o (r r 2 ) 



(m + l)R e 

p | Dm+l pm+1 

I A e,max 

(m + 1) f 1 - ft 



o (r- 2 ) 



15) 



+ o (r- 2 ) 



f m (Rt)-RT +1 



(m + 1) ( 1 - ft 



i + o (r- 2 ) 



i6) 



and for the derivatives 

dfle (m + l)R™ +1 { !-/?,) ! /,„(/?,)-/?,'" 1 



eiiL 



2(m + l)r?ftf 



+ o (r- 4 ) 



(87) 



dAir (m + l)^ 1 (l " Re) - [URt) - 



dR P 



2(m + l)r t 2 1-ft 



o(r- 4 ) . (8? 



The Doppler factor is given by 



2T 



2(m + l)T t R 



(m+2)/2 



1-it 



l + r 2 (^ e +^) 2 (7^ + 1)^+1(1-^)+ f m {R t ) ~ R? +1 



(89) 



where we have used eqs. (|85|) and (1561) as well as T 2 = TfR e m 3> 1 and 9 e + 9 r <C 1. Finally, 
f ~ 1 — -R e , and thus 

-(l+a) 



5 Q+2 d/i e ^5 Q+2 d/i £ 



2(2r 4 ) c 



f 2 dR 1 f 2 ATi (-\ 6 



r| dft e (1 - # e )i^ 



1 



fm{Rt) - RT +1 

(m+ 1)^+1(1 -ft e ) 



2(2r t ) a (m + i) 1+a ft; 



l+a p a + m ( 2 +a)/2 



(l-^Re)" 



m + l)^+i 1 - Re + / m (ft t ) - ft e m+1 



l+a 



(90) 
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We note that the expressions above for case 3 are identical to those for case 1 - Eq. (j73 



\6 r — 6 T 



-1/2 



at fl^max, dfli/ d,Ei 



While drii/ deidVL r diverg eS as \R e — R e [yr max) I 
J dQ r (drii/deidQ r ) ~ 2tt J 6 r dO r (drii/dEidQ r ) remains finite (i.e. both the energy density and 
the energy flux of the radiation field remain finite). This has been noticed in the context of 



the diverging surface brightness of the afterg 
comes from an in finitely thin she l l dSari Il998 



ow image at its outer edge, when the emission 



Granot fc Loeb 120011). In that context, it has 



also been shown (jWaxmanl 119971 ; iGranot. Piran fc Saril Il999a| jbl; iGranot fc Loeb! 120011 ) that 
when the emission comes from a shell of finite width, the surface brightness (i.e. the specific 
intensity I £ ) does not diverge. 



4.4. Putting it all together 

Analytical expressions for our model have now been fully derived, and are reported for 
convenience here. The scaled spectral flux density, Eq. ( fl"3l) is rewritten as: 



26— rna -i , „ 

i+a r ' AR T 



F £ (T) (T\^) ( m + 1 » ^!_ ma / 2 

d y [ —— _ m _i 1 V ma/2 exp 



F £ ,o \T J J y . \m + y 



Ro ' To 



(91) 

where e t = (1 + z)e, y min = min[l, R Q /R L (T)] and y mSuX = min[l, (Rq + AR)/R L (T)], while 
the flux normalization is given by 

(1 + z) _ 2 a L e 1 ~ a (l + zf~ a 
F Efl - 2T L {1+z)£/2To (R ) - ^ , (92) 

F = s^F^o = (sF £t0 ) \ E=1 = A (^) (93) 



ndi V 2 



i \ 2-a 
1+2 



- 7.6 x lO" I ^ j L ^ 2 dl% erg cm" 2 s" 1 , 

where d^ = lO 28 ^^ cm ; an d ma Y be used in order to infer the value of L from the observed 
flux level. The optical depth in the integrand above is: 

r 77 (e t , e tt0 , Rt, ) = n i?J~ 6 ~ ma/2 .F(x) , (94) 

where Rq = y~ l (T /Tq)' 1 /^ 1 ^ and x = (y-( m+1 ) - l)/(m + 1). The function T is the 
following double integral: 



F{x) = dR t I dR e X(Rt, R e ) + / dR t dR e X(R t , R c ) . (95) 

J l JRo/Rt JR2 JRo/Rt 



-28- 



The two integrals above correspond to cases 1 and 3 respectively, as discussed in § 14.31 
When x > 1, the test photon lags behind the shell and f m {Rt) > 1, with f m (Rt) = [1 + 
x(m + 1)(1 — Rf 1 )]/ R r ^ l+1 ; R 2 is then defined by the implicit equation f m (R 2 ) = 1, and 
Re t 2 = min[(i?o + AR)/R t , 1]. For m = 0, R 2 = max(l,x), while for m = 1, R% = 
max[l, + 8x — l)/2]. The test photon eventually overtakes the shell at radius R 2 , and 
will travel ahead of the shell ever after, which corresponds to the second integral, where 
-Re )3 = min[(R\, + AR) / R\, f m {Rt) 1 ^ mJrl \ When x < 1, only the second integral contributes, 
as the photon is emitted on the shell and immediately travels ahead of it. Thus, R 2 — 1 so 
that the first integral vanishes. Note that, for all practical purposes, Eq. fl95l) can be cast 
into a single integral: 



dR, 



dR e I[Rt, R e 



(96) 



Ro/Rt 



where R, 



R 



i? e> 2 when x > 1 and f m (Rt) > 1, and i? e ,M = -R e ,3 in all other cases, i.e. 



e,M 



mm[(R + AR)/R t , f m (R t 



T(Rt, R e 



5 3 d\i e 



R? \r 2 dR, 



1 V(™+1) 



1]. Finally, the integrand is equal to: 



c- V(C+-C)(C-C-) h x °t 



L\,2 



X 2 /et(S 



(Re)- (97) 



Specializing to L' E ,(R) = L'e/ x (R e R t /R f, Eq. ® becomes: 



i)-2+mo/2 



dR 



Ro/Rt 



dR 



where the integrands are further expressed as: 

P +a dfi e ~ b [2(m+l)] l+a R b e +a+ ^ {2+a \ 



-=--tl„ 
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(98) 
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m + l)^ 1 ( 1 - R e ) + / m (i2t) - i? e m+1 

4(r t r )(r t e t ) 



1+a ' 
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(m+l)(l-i? e ) 
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H a (()= 2 F 1 (-a,0.5;l;-() , (99) 

(100) 
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This concludes the set of general equations that have been obtained. For reference, the 
hypergeometric expressions for a — 1, 2, 3 respectively read 



Ci#i(C) = C- i + 



c 



1 



[(T t e r ) 2 + (T t e t ) 2 ] , 

(v t e t ) A + (r t e r y_ + (r^) 2 (rA) 2 



16 

c 3 ^3(o = e(i + ^+^ 2 +^c 3 

= ^ [(W 6 + 9(r^ r ) 4 (rA) 2 + 9(r t 9 r ) 2 (r t 9 t ) 4 + (r^) 6 ] 



For our fiducial case, a = 2, we also explicitly write the relevant expressions : 

-4 

m-6+i ^ _ n /i no / 1 \ 



r (e t , R tfi ) = r* e t i? p - m - 6+1 , t+ = 0.402 
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m + l)(l-i? e ) 16 (m+l) 2 (l-i? e ) 2 
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(104) 
(105) 
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5. Analytic Scalings of the Flux and Optical Depth 

Before showing our results for the lightcurves and spectra, it is useful to first analyti- 
cally derive some of the relevant scaling laws (from the equations obtained in the preceding 
sections), and discuss the qualitative behavior of the system in different regimes. It is con- 
venient to define a normalized time T = (T/T ) — 1, which is zero when the first photon 
from R t = R Q and 9 t = reaches the observer, and is ~ 1 about a dynamical time later, 
when the system starts to approach a quasi-steady state. It is also useful to define the time 
T f = T (l + AR/R ) m+1 , where R L (T f ) =R + AR, when the lack of emission from outside 
the outer edge of the emitting region (R > R + AR) starts being noticed by the observer, 
and the corresponding normalized time 

( (m + 1) AR/Ro < 1 (AR < R ) , 
- 1 « I (108) 
{ (AR/R ) m+1 > 1 (AR > R ) . 



T f ( AR\ 

T ' = i- 1 = { 1 + -^) 
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Note that at T < Tf the outer boundary of the emission region does not affect either 
the emission, since the outer edge of the EATS-I is still fully within the emission region, 
or the opacity of the emitted photons, since the maximal radius of the EATS-II (i£ e ,max) 
at all points along the trajectory of any photon is always smaller than that of the EATS-I 
[Rl(T)}: i? e)max (T, y, R t ) < Rl(T). As shown in Fig. [H the two radii become nearly equal for 
R t Rl(T). In fact, for R t 3> Rl{T), not only does R e)mBX (T,y, R t ) approach Rl(T), but 
the EATS-II approaches the EATS-I (the two must become identical when the test photon 
reaches the observer, which corresponds to R t — > oo for a distant observer, at "infinity"). 
This immediately implies that for T < Tf the observed flux and the opacity along the 
trajectory of all photons (which reach the observer at time T) are independent of AR. Thus, 
in order to calculate the light curves for a family of model parameters that differ only in their 
AR values, it is sufficient to calculate the observed flux and opacity for AR — ► oo and use 
them for T < Tf, and do the full calculation for each specific value of AR only for T > Tf. 

The temporal scaling of the unattenuated flux, at sufficiently low photon energies e, 
can be understood as follows. For 1 > f < f f , y min = R /R L {T) = (T/T y l ^ m+1 ^ = (1 + 
y^-i/(m+i) anc [ y max — x j go that Ay = y max —y m i n ~ T / (m+1) while y ~ 1 and the integrand 
in Eq. (1131) is also ~ 1, implying that F e oc T. For Tf < T <C 1 the emission is from R m Rq 
from angles 9 which satisfy {To9) 2 ~ Tf (m + 1) <C 1, while the Doppler factor at this stage is 
almost constant, 5 w 2T /[1 + (T 9) 2 ] w 2T /[1 + f/(m+ 1)] w 2r , and therefore the flux is 
approximately constant in time, F £ oc 5~ l ~ a oc T°. For Tf <^ 1 <^ T, 5 ~ 2r /(r 6') 2 oc T _1 
and F e oc T _1_a . For 1 < f < f f , y mhl = R /R L (T) = (1+f )~V(™+i) ~ f-i/(m+i) < i an d 
Umax = 1, so that the integral over y in Eq. ffTBT) approaches a constant (corresponding to its 
value for J Q dy), and F e oc j'( 26 - ma )/[ 2 ( m + 1 )] , Finally, for T Tf, the emission is dominated 
by R ~ Ro + AR and an gles 6tp S> 1/r^n (i.e. x 3> 1 ) and we obtain the familiar result for 
"high latitude" emission flKumar fe Panaitescull2000f ). F e oc T~ 1_a . Altogether, 



-^e<ei(T) OC < 



T (1 » T < T f ) , 

T° (27 < f < 1) , 

j.(26-ma)/[2(m+l)] (1 < T < Ty) , 



-1-Q 



(109) 



(T > max[l,7>]) 



Now we move on to discuss the opacity effects in some detail. As can be seen in Fig. 
at a given emission radius the optical depth is smallest for small emission angles (i.e. small 
values of x). There is a local maximum near x = 7^0^.0 ~ 1 since for such emission angles the 
photon is emitted almost parallel to the shell in the comoving frame, and a relatively large 
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part of its trajectory (also in the lab frame) is close to the emitting shell, which enhances the 
optical depth. For a given normalized emission angle, x 1 ^ 2 = jt,o8t,o, the normalized optical 
depth increases with emission radius, as can be seen in Fig. where the increase is largest 
for small emission angles. The optical depth generally increases with AR/R when all other 
model parameters are held fixed, due to the larger range of emission radii which enhances the 
photon field that can potentially interact with test photons. However, as pointed out above, 
for T <Tf the optical depth in this case is independent of AR/R . This is demonstrated in 
the lower panel of Fig. [HI where it can be seen that in practice a noticeable increase in the 
optical depth due to the increase in AR/Rq does not occur immediately after Tf but takes 
some time to come into effect. This is since for < T — Tf <C Tf the added contribution 
to the opacity from R > Rq + AR for the smaller AR is very small, since the additional 
photons can interact with the test photon only at very large radii (R t 3> Rl(Tj) = Ro + AR) 
where the intensity of the photon field is very small, and at very small angles between the 
directions of the photons which are very unfavorable for interaction. 

One would also like to define S\ as the photon energy at which the optical depth becomes 
unity: r 77 (£i) = 1. However, this definition gives a different value along the trajectories of 
different (test) photons, making it hard to define a unique value for £i(T), since its value 
varies along the EATS-I (see Fig. [7]). For 1 <f < Tf, F(x) becomes independent of T and 
depends only on x (see upper panel of Fig. [7]). Most of the contributions to the observed 
flux come from x < 1, since for x ^> 1 the radiation is strongly beamed away from the 
observer. The upper panel of Fig. [7J shows that Tix < 1) varies over a factor of ~ 50 
for 1 <C T < Tf, and therefore in this regime it still makes some sense to define a single 
typical value of £i(T), and derive its scaling. It is good to keep in mind, however, both the 
spectral transitions around ei(T) in the instantaneous spectrum and around £i(Tf) in the 
time integrated spectrum, as well as the transition in the light curve when £\{T) sweeps past 
the observed photon energy e, are all expected to be somewhat smoothed due to this relatively 
large range of opacity values across the (unresolved) observed image of the GRB projected 
in the sky. According to Eq. ([38]), r 77 oc r oc e t a ~ l R\~ h ~ mal2 oc e t a ' x f -(i-&-W2)/(m+i) ? 
and therefore 6i(l«T< Tf) oc r(i-&-"W2)/[(m+i)(a-i)] _ 

The lower panel of Fig. [7J shows T[x) as a function of Y = (y — y m m) / (Umax — Z/min) ~ 
(x max — x)/x max for several values of 1 ^> T < Tf, along the equal arrival time surface of 
photons to the observer (EATS-I). In this limit R tt o ~ Rq, y ~ l — x and x max ~ T/(m+l) ~ 
[r -R_L i max(T')/-Ro] 2 where R± t „iax(T) is the radius of the GRB observed image, projected in 
the sky, at a normalized observed time T. As is shown analytically in Appendix [D] and is 
apparent in the lower panel of Fig. [7J in this limit 



T{x = 0, 1 > T < T f ) oc T 



(110) 
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The lower panel of Fig. [7] also shows that 

( Y Y> y,(T) , 
^(s, 1 > T < T f ) w J^(x = 0, 1 > T < 2» x ^ (111) 

ln(f) F<Y»(T), 

where Y* ps (x max — x*)/x max oc T a is the value of K where the scaling of J-{x) changes from 
oc Y° to oc Y . The corresponding value of x is x*, and ~ f(i mM ) oc 

The scaling of T{x) with Y" for 1 ^> T < Tf can be understood as follows. For a given 
emission radius i? i , the dependence of the optical depth and J-"(x) on the emission angle 
t fl is very weak for x = (r t)O t)O ) 2 ^ 1, and becomes significant only for a; > 1 (see Fig. [5]). 
Therefore, as x starts increasing from x = at the line of sight, along the EATS-I, !F(x) 
initially varies following its dominant radial dependence. The latter may be derived from 
that along the line of sight where T(x) oc T oc (R tt0 /R Q ) — 1, where for a general value of x < 
i max < 1 we have (Rt t0 /Ro) — 1 ~ x max — x ocY and therefore J-{x) oc Y. As x approaches 
i max , Y approaches zero, until eventually the optical depth becomes dominated by the small 
angular dependence on 6 t $ at a fixed emission radius R tt o, and F{x) approaches a constant 
value of f(x mm ) which corresponds to R t) o = Ro and x = (T 9 tt o) 2 = x max = T/(m + 1). 
We find numerically that f(s max ) oc T a+1 oc = (r 6 l 4j o) 2( ' a+1 - ) . This may be understood 
as follows, starting from the expression for the optical depth in Eq. fl3Tl) . In this regime 
-R e ,max = /m (m+1) and for R t = R we have 



-Re, max ^ 



-i l/(m-l) / „ \ 

l + x(m + l)(l-i?7 1 ) - 1 ^x fl -i^M < 1 • (112) 



This means that the contribution to the local photon field at each point along the trajectory 
of the test photon is always from a very narrow range of radii near R . This implies that 
r~~ 2 \d/j, r /dfj, r \ which appears in Eq. fl3T|) remains approximately constant, since the geometry 
of the problem implies Rq6 c ^ r6 r so that r~ 2 \dfx r /dfx r \ ~ R$ 2 = const. At any given 
point along the test photon trajectory r < 9 t < 9t,o, simply because in this regime ^ max 
is obtained where EATS-II is truncated at Rq, which must correspond to 9 r = 8t for a 
test photons that is emitted at R t $ = Rq (the test photon is always on its own EATS-I 
and EATS-II, by definition). This implies that 5 ~ 2r = const since T(6 e + 9 r ) < 2T 6 t < 
2r #to = 2a; 1 / 2 <^ 1. Furthermore, L' ,(R e ) is approximately constant since R e ~ Ro = const. 

i 

Since 9 t /9 t ,Q ~ Rtp/Rt = Rt, the effective solid angle that contributes to interaction at R t is 
~ d\ oc R~^ 2 and there is also a factor of 1 — \i ti ~ 9 2 oc Rj 2 in the integrand of Eq. (l3Tj) . most 
of the contribution to the total optical depth is from R t ~ R tt0 (i.e. R t < 2). Therefore, the 
integration over the solid angle effectively introduces a factor of ~ 9 2 , while the factor of 
1 — [L t i in the integrand introduces a similar factor, together giving a factor of ~ 9f Q . The 
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integration over energy Ei together with the threshold e t Si > 2/(1 — /j, ti ) for pair production 
give L' £ , oc e\~~ a oc £? -1 (l — ^u) 01 ' 1 oc Ofja ■ Altogether, with the previous factor of 8f , 



the optical depth in this regime scales as 6. 



2(a+l) 
t,0 



Thus, for fixed values of r* and e, r 77 first becomes larger than unity at the center of 
the image (i < 1 and Y ~ 1) at Tij ~ a 2 2a_1 /r^e Q_1 . From this time on the central part of 
the image is opaque, at x < x\ which corresponds to Y\ « Tu/T, so that photons of energy 
£ can escape mainly from a thin ring in the outer part of the image, that corresponds to 
x\ < x < x max and occupies a fraction Y\ ~ Tu/T oc T _1 of the image area (since that area 
is linear in x). Thus, the observed flux is suppressed by a similar factor and turns from oc T 



at T < Tu to oc T° at T > Tu. Eventually, at a later time T\f ~ 7\/ oc e 
when Xi 



-(a-l)/(a+l) 



3? * 



the whole image becomes opaque, i.e. r 77 > 1 for all < x < x 



and 

the observed flux starts to drop exponentially with time. This behavior can be seen e.g. in 
Fig. [8l In summary, 



F £>>£l ,{f <l)~F e<ei .(f = l)x{ 



T 



■ it 



f^exph^/f^) 



a+l] 



f u {e) < f < f lf {e) , (113) 
f > f lf (s) . 



Similarly, for 1 ^> T < Tf it is natural to define £u(T) and £i/(T) as the two photon 
energies above which the center and outer edge of the observed image, respectively, become 
optically thick to pair production: by definition, Tuj[su,f(T)] = T. This implies that 
su ~ (a2 2a - 1 /nf) 1 /( a -V oc r-V(«-i), and since f x/ ~ T 1 Y ( ' a+1) , we have e u /e lf ~ f "Z^ 1 ) 
and £1/ oc J , ~( a+1 )'( a ~ 1 ). Eq. (jl 13j) determines the instantaneous spectrum in this regime, 



F e>ei ,(T < 1) ~ TF e=1<£l XT =l)x{ 



-(a-l) 



a-1 -2(a-l) 
fc li fc 



£ < eii(T) , 
£ H (T ) < e < e lf (f ) , 



a-l -2(a-l) 



exp[-(e/ei/) 



a-ll 



£ > £i/(T) . 

(114) 

At T ~ 1 the opacity becomes more uniform across the image, Tu ~ Tf ~ 1, and Eu ~ 

Eif ~ £i(T = 1) ~ £i*. 

For e ^> £i(T = 1), the time integrated flux f £ = J dTF £ (T) is approximately given by 
~ T F £<£l (f = l)f u f lf where f^f^ oc 7\ ( f +2)/(q+1) oc e -(a-i)(«+2)/(«»+i) > since f u oc e 1 -". 
Therefore, the spectral slope of the time integrated spectrum, f £ , for impulsive sources 
(AR < Rq and Tf < 1) where the total time integrated flux is comparable to that from the 
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rising phase, steepens by Aa = (a — l)(ot + 2)/ (a + 1) above ei(T/), 

-(a-l)(2a+3)/(a+l) [ £ > 



/ e (Ai2 ~ i? ) oc 



(115) 



This can be seen in Fig. [121 For a quasi-steady source (Ai? ^> _R and Tf ^> 1), a similar 
time integrated spectrum is obtained only if the flux at 1 < T < Tf decays faster than T -1 , 
i.e. if m(a — 2) > 2(6 + 1) [see Eq. (11091) ]. so that f £ is dominated by contributions near 
T ~ 1. For a slower decay or a rising flux at 1 < T < Tf, f e is dominated by contributions 
from T ~ Tf and there is an exponential cutoff above ei(Ty), while the power law high 
energy tail from the rising phase is encountered only after a significant (exponential in e) 
flux drop. For extremely impulsive sources, where Tf <C 1 (i.e. AR <C R ), there is also an 
intermediate power law segment in the time integrated spectrum: 



-(«-i) 



< eu{f f )] 



f £ (AR « Rq) oc { e- 2 ^ [e u {ff) < e < e lf (T f )} 

.-(a-l)(2a+3)/(a+l) [ £ > £l/ (f f )] . 
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6. Results: Semi-Analytic Light Curves and Spectra 

Figures [8] - [11] show light curves and spectra for the semi-analytic model developed in 
the preceding sections. We use fiducial parameter values of m = b = 0, AR/R Q — t* = 1, 
and a = 2, which are relevant for the prompt gamma-ray emission in GRBs, and vary one 
parameter at a time in order to see the effect of each model parameter more clearly. When 
varying m and b (Figs. [9] and [10l respectively) we use AR/R = 100 in order to have a 
large enough range of emission radii so that the radial dependence of the Lorentz factor and 
of the co-moving spectral emissivity would have a significant effect on the light curves (for 
AR/ Rq <^ 1 the values of m and b hardly affect the light curves). Figure [121 shows the time 
integrated spectra for several values of AR/R , where each panel is for a different set of 
values for the three parameters (a, m, b). In order to ease the reading, Table [2] summarizes 
the various sets of parameters and the corresponding figures. 

Figure El shows the light curves for fixed values m = b = 0, r* = 1, and a = 2, while 
the various panels correspond to different values of AR/Rq (of 0.01, 1, and 100, from top 
to bottom). At the lowest photon energies, well below eu (which for the parameter values 
used here is ~ 10 2 ), opacity to pair production never becomes very significant, and the light 
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curves follow the behavior described in Eq. (11091) which is discussed in the preceding section. 
In this regime the lightcurves are self similar in the sense that e a ~ 1 F £ is independent of e 
below £\{T). The different behavior for Tf <C 1 and Tf ^> 1 (where Tf = AR/Rq for m = 0) 
that appears in Eq. (11091) can clearly be seen by comparing the upper and lower panels of 
Fig. [HJ For e ^> Eu, on the other hand, opacity to pair production has a major effect on the 
light curves. In this regime the light curves at T <C 1 follow Eq. f 1 1 1 3 [) . showing a pronounced 
constant flux plateau between Tn oc £ 1_a , when the center of the image becomes optically 
thick to pair production, and T\f ~ T\[ ^ a+1 ^ , when the entire image becomes opaque, followed 
by an exponential flux decay. At 1 < T < Tf the opacity does not vary drastically across the 
image and may be described by a single value of ei(l <C T < Tf) oc j ] ( 1 ~ b - ma / 2 )/[( m + 1 )( a - 1 )} . 
For the parameter values used in Fig. [HJ E\ increases (linearly) with T in this range, and 
therefore the opacity at a given e decreases with time, causing the observed flux to increase 
with time until E\ sweeps across e or until Tf is reached (whichever comes first). At T > Tf 
the situation is reversed, as the observed emission comes from large angles relative to the 
line of sight ( "high-latitude" emission) and E\ decreases with time. 

We now turn to the photon energy spectrum. The instantaneous spectra at T <C 1 follow 
the behavior described in Eq. (11141) . At very early times the exponential part starts only 
at very high photon energies, making it very hard to detect. When T ~ 1 the intermediate 
power-law segment disappears as en and e\f become nearly equal (note that the low energy 
part of the curves appears flat in the figures since we show e a ~ 1 F £ which is independent of e 
below Si). The time integrated spectrum varies with the value of AR/Rq. For AR/Rq <C 1 
it consists of three power-law segments, as described in Eq. (jll6p . As AR/Rq increases, the 
central power-law segment, at Eu(Tf) < e < Eif(Tf), shrinks as £n(Tf) and E\f{Tf) approach 
each other, until it disappears for AR/Rq ~ 1 where Tf ~ 1 and EuiTf) ~ E\f{Tf) ~ Eu. 
For AR/Rq ~ 1 the time integrated spectrum is described by Eq. (I115p . and consists of 
two power-law segments. As AR/Rq increases above unity the time integrated spectrum 
develops an exponential high-energy cutoff, while the power-law tail at high energies becomes 
increasingly suppressed. This occurs since if the flux at 1 < T < Tf does not drop faster than 
T _1 , which corresponds to m(a — 2) < 2{b + 1) (see Eq. [109] ) as is indeed the case for the 
parameter values used in Fig. [SJ then the time integrated flux is dominated by contributions 
from T ~ Tf ^> 1 and reflects the exponential cutoff of the instantaneous spectrum at 
that time, which dominates over the high-energy power-law component that arises from the 
superposition of the instantaneous spectra from T < 1. 

Figures M and [TD] demonstrate the effects of the two parameters m and b. As discussed 
above, a large value for AR/R (100) was chosen so that the radial dependence of the Lorentz 
factor (r 2 oc R~ m ) and of the co-moving spectral luminosity [L' £ , oc R b (E , ) 1 ~ a ] would have 
a large effect on the light curves. For AR/Rq 1 the values of m and b hardly affect the 
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light curves (since the emission takes place over a very small range of radii in which both T 
and L' e , hardly vary). Figure [9] also demonstrates the dependence of Tf on m, where in the 
limit of AR/R > 1, T f » (AR/R ) m+1 (see Eq. [TUS] ). 

As can be seen in Fig. [9] the power-law component of the time integrated spectrum 
is largely independent of m, since it originates from the superposition of the instantaneous 
spectra at T < 1, which are sampling a small range of emission radii. The lower energy 
component, however, from the contribution of the emission at times 1 < T < Tf, is sensitive 
to the value of m, since it sample a large range of emission radii. For m — 0, F £ (l < T < Tf) 
is constant in time (for the values of the other parameters that are used in Fig. [9]), while 
£i(l < T < Tf) oc T, and both effects combine to produce a very pronounced high-energy 
exponential cutoff. For m — 1, F e (l < T < Tf) oc T -1 / 2 while £i(l < T < Tf) is constant 
in time, which results in a somewhat less pronounced, though still fairly large high-energy 
exponential cutoff in the time integrated spectrum. For m = 2, F £ (l < T < Tf) oc T~ 2 / 3 
while £i(l < T < Tf) oc T^ 1 / 3 , so that the time integrated spectrum in the range £i{Tf) < 
e < Sx(T = 1) ~ Su is dominated by the contributions near the time T(e) when £i(Ti) = e. 
This results in a spectral slope of ef £ oc e~~ l in the lower panel of Fig. [9j 

More generally, F £<£l{T) (l < f < f f ) ~ e i- a f{2b-m a )/[2{ m +i)\ while £l (\ < f < f f ) ~ 
j , ( 1 - ( '- ma / 2 )/[( m + 1 )( a - 1 )] ) go that when the flux is dominated by the contribution from 
T ~ Ti(e), then the spectral slope of the time integrated spectrum is given by 

dloge^fs = (a-l)[m(2-a) + 2(6 + 1)] 
rfloge 2(1 — b) — ma 

This may be relevant if Ei(l < T < Tf) decreases with T, in which case this spectral slope 
is valid in the range £\{Tf) < e < e^. It may also be relevant if £i(l < T < Tf) decreases 
with T, as discussed below. 

In Fig. [TDJ, the upper panel is identical to the upper panel of Fig. [9] and the lower panel 
of Fig. E In the middle panel F £ (l < T < 7>) oc f" 1 and ei(l < f < f f ) cx f 2 , while 
in the bottom panel F £ (l < T < Tf) oc T~ 2 and 6i(l < T < Tf) oc T 3 . Both cases result 
in a very pronounced exponential cutoff at very high photon energies (which may be hard 
to detect), but show a shallow spectral slope up to this exponential cutoff (which may be 
easier to detect). This again results in the spectral slope given by Eq. (11171) . However, in 
this case Ei(l < T < Tf) increases with T, and therefore this spectral slope occurs in the 
range < e < 6i(Tf). This is valid, however, only if indeed the time integrated flux in 
this spectral range is dominated by the contribution from near Ti(e). This is not valid in 
the upper panel of Fig. [10] (where it is dominated by the contribution from T ~ Tf), and is 
only marginally valid in the middle panel (where the contributions from all the times in the 
range Ti(e) < T <Tf are comparable). In the lower panel of Fig. [10] the flux in this spectral 
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range is indeed dominated by the contribution from T ~ T\(e), which results in a spectral 
slope of ef £ oc e~ 1 ^ in this range. 

In Fig. [HI which shows the effect of varying a, the top panel corresponds to a = 1, for 
which both the flux and the optical depth become independent of e. As a result, we present 
for this case light curves for different values of t*, and the corresponding integrated spectra 
(which all have a flat f £ and vary only in their normalization). For a = 2 (middle panel) and 
a = 3 (bottom panel), one can verify that the power laws on the middle and bottom right 
panels have an index of approximately 4/3 and 5/2 respectively, as expected from Eq. (11151) 
after rescaling by e a ~ x . 

Finally, Fig. [12] illustrates the behavior of the time integrated spectra, as discussed at 
the end of § [5J All the curves show a high energy power law tail with an index of about 
4/3, as expected from Eq. (11151) . Moreover, given the rescaling by T (m = 0) (which is 
independent of m) in Fig. [121 it is easier to see that the time integrated spectra become 
independent of b and m for AR/Rq <C 1 (since in that limit, the same holds for the light 
curves and instantaneous spectra). As discussed in the paragraph following Eq. (11151) . the 
exponential cutoff is suppressed when m(a — 2) > 2(6+ 1), as is the case on the middle panel 
only. In such a case, the time integrated spectra are dominated by contributions near T ~ 1, 
and the effect of AR/Rq becomes negligible for AR/Rq 3> 1, which explains the asymptotic 
behavior of the spectra with increasing AR/Rq. Finally, for very impulsive sources, the 
intermediate power law segment in Eq. (11161) can be discerned, albeit with difficulty 

7. Discussion 

We have explored in great detail a model for the temporal and spatial dependence of 
the opacity to pair production (77 — > e + e~) in impulsive relativistic sources. Our simple, 
yet rich, model features a thin spherical shell expanding ultra-relativistically and emitting 
isotropically in its own rest frame within a finite range of radii. Our two main results are the 
follwoing. First, while the instantaneous spectrum (which is typically very hard to measure 
due to poor photon statistics) has an exponential cutoff at high photon energies, the time 
integrated spectrum over the duration of a flare or spike in the light curve (which is easier 
to measure) has a power-law high-energy tail. Second, photons above this spectral break in 
the time integrated spectrum arrive mainly near the onset of the flare or spike in the light 
curve. 

These two features provide a unique detectable signature of opacity to pair production, 
making it easier to identify observationally. Furthermore, these features are expected to 
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be fairly robust, even if the exact details (such as the exact change in the spectral slope 
across the break, Act, or the exact shape of the light curve at high photon energies above 
the spectral break) may depend on the details of the model (such as the exact geometry, 
which is assumed to be spherical in our mode|§). The reason behind these features is that 
in impulsive sources the photon field starts from zero (or more realistically a non-zero value, 
which is still much lower than that near the peak of the flare or spike in the light curve) and 
builds- up with time, so that the optical depth to pair production, r 77 , increases with time, 
and high energy photons can escape mainly at early times while r 77 is still below unity. 

A source is considered impulsive for our purposes if the photon field in the source and 
its vicinity changes considerably within the source light crossing time. In this limit the time 
dependence of the photon field and the resulting opacity to pair production is important. 
This can naturally occur in relativistic sources, but is hard to produce in non-relativistic 
sources (since it requires a relativistic signal in order to turn the emission on or off on 
a time scale of the order of the light crossing time of the source). In the opposite limit, 
where the photon field hardly varies within the source light crossing time, the photon field 
may be approximated as constant in time along the trajectory of the photons, and can 
be evaluated at the time of emission (this is considered a "quasi-steady" state). In our 
model, the photon field approaches a quasi-steady stat^ within a few light crossing times 
of the emitting region (Tf > T > a few). If the source is active for much longer times 
(Tf > 1 <o> AR/R ^> 1), the time integrated spectrum will usually be dominated by this 
late time quasi-steady emission, and an exponential high-energy cutoff develops, while the 
high-energy power-law tail becomes increasingly suppressed. For this reason, in order for 
the source to be impulsive, the duration of the emission should be at most comparable to 
the light crossing time of the source (Tf < 1 AR/Ro < 1). 

We have considered a single, isolated emission episode which corresponds to a single flare 
or spike in the observed light curve. Furthermore, we have assumed no background photon 
field at the time when the emission turns on. These are obviously idealized assumptions 
and it is worth considering, at least qualitatively at this stage, the modifications that may 
occur when these ideal conditions are not satisfied. For the prompt emission or X-ray flares in 
GRBs the background quasi-steady photon field is expected to be very low and not contribute 
significantly to r 77 . This is probably also true for Blazars of BL Lac type. In other types of 
Blazars (quasar hosted Blazars), however, the external photon field, mainly due to emission 



5 In AGN, e.g., a cylindrical geometry may be more appropriate. We intend to study such a cylindrical 
geometry in a future work. 

6 Here, by quasi-steady state, we mean that neglecting the time dependence of the photon field would at 
most change the results by a factor of order unity, but not qualitatively. 



-39- 



from the accretion disk and its scattered photons from the clouds in the broad line region, 
i s expected to provide the dom inant contribution to r 77 in the source, even during flares 



( ISikora. Begelman Sz Reeall994j ). In this case, our model would not be applicable, since the 
external radiation field would prevent the escape of high-energy photons near the onset of 
the spike, resulting in r 77 that is largely independent of time. 

Regarding the assumption that the flare/spike is isolated, in many cases there are series 
of flares, so that except for the first flare in the series for which our assumption should hold 
very well, for consecutive flares the high energy photons could in principle pair produce 
with photons emitted in previous flares. This will be highly suppressed if the time from 
the end of the previous flare is much larger than its duration. Even if these two times 
are comparable, pair production is still significantly suppressed since the relevant photons 
can meet only further away from the source, where the photon density is smaller and the 
angle between the directions of the photons (in the lab frame) is smaller. Such a proximity 
in time to a previous flare/spike will still increase r 77 to some degree, but this will affect 
mainly the highest energy photons, with energies well above the spectral break in the time 
integrated spectrum over the duration of the flare/spike, which are relatively hard to detect 
due to the smaller number of photons at such high energies. Therefore, in practical terms, 
the differences from our idealized model are not expected to be very large. This may even 
make it meaningful to integrate the spectrum over many spikes/flares in order to increase 
the photon statistics, in cases where the number of photons detected in individual spikes is 
not large enough to enable a good spectral analysis. 

Other sources of opacity, such as scattering of photons on the pairs that are produced, 
are also possible. The latter, however, is expected to build-up in time on a comparable 
time scale to that of the opacity to pair production that we study here. Therefore, it is 
not expected to have a significant impact on our main conclusions. Opacity for scattering 
on the electrons associated with the baryons or with preexisting pairs within the outflow 
is also possible. However, it will not greatly vary within a single dynamical time, and 
should also affect lower energy photons (where the larger number of photons enables a 
better spectral analysis). Furthermore, it is suppressed at high photon energies due to the 
reduction in the cross section in the Klein-Nishina regime. Moreover, we find a rather unique 
combined spectral and temporal signature for the opacity to pair production, which could 
help distinguish between it and other sources of opacity. 

In GRB afterglows the opacity to pair production is typically very low and therefore 
not expected to be detectable in the GLAST energy range. During the afterglow, after 
about one day, £0,52 ~ 10~ 8 — 1CT 7 , R ~ 10 17 cm corresponding to -Ro,i3 ~ 10 4 , and 
r ~ 10 corresponding to r ,2 ~ 0.1. According to Eq. f)12ip . this implies a huge value of 
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£1* (ei*m e c 2 ~ 10 15 — 10 16 eV for a ~ 2). In practice the opacity would be even lower than 
this, since the typical energy of the photons that would pair produce with such high energy 
photons would be ~ r 2 /si* corresponding to v ~ 10 12 Hz, which is well below the assumed 
power law segment of the spectrum (so that the number density of these low energy photons 
would in practice be much lower than the default value according to our assumption of a 
simple single power law spectrum). A possible exception to the very low r 77 during the 
afterglow may be the very early afterglow in a stellar wind environment, near T^ cc which 
is of the order of seconds in this case. Typical parameters values there are -Ro,i3 ~ 100' 
r ,2 ~ 1, and Lo,52 ~ 0.1, which might give eum e c 2 as low as ~ 100 GeV. Such values could 
be detected by GLAST, albeit with difficulty 

For the internal shocks model, the GRB prompt emission occurs at a much smaller 
radius compared to the afterglow, -Rqrb "C -Rdec- Furthermore, the luminosity of the prompt 
GRB emission is much larger than that of the afterglow, and the Lorentz factor is higher. 
Therefore, despite the higher Lorentz factor during the prompt GRB, r 77 is still much larger 
than during the afterglow. In models where the prompt emission occurs near the deceleration 
radius, -Rgrb ~ -Rdcc, the values of r 77 in the very early afterglow and in the prompt emission 
are comparable (perhaps somewhat smaller in the early afterglow due to a smaller radiative 
efficiency), but r 77 is typically very low for both emission components (i.e. the effects of 
opacity to pair production are not expected in the GLAST energy range). 

Once the spectral break in the time integrated spectrum over the duration of a flare 
or spike in the light curve is observed in the data, it can be used to constrain the values of 
the physical parameters of the source, namely Tl a R . A fit of our model predictions to the 
data can in principle determine the values of all the model parameters: a, m, b, AR/Rq, 
and F , which in turn determine L (from F , Eq. [93]) and T^Rq (from r*, Eq. |35j). 
In practice, however, the limited photon statistics may render such a direct fit with such a 
large number of free parameters impractical. One way to overcome this problem is to fix the 
values of some of the model parameters, e.g. m = b = and even AR/Rq = 1, if necessary. 

A less accurate but less computationally demanding alternative is to fit the time in- 
tegrated spectrum (over a flare or spike in the light curve) to a parameterized function 
featuring a smooth transition between two power laws, 



where n and fo determine the sharpness of the spectral break at Eu and its flux normalization, 
respectively, while / e < £l „ oc e 1 ^ and / e > ei , oc s 1 ^^^. Such a fit can determine both the 
photon index, a, and the photon energy ~ S\{T = 1) of the spectral break in the 




l/n 



(118) 
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time integrated spectrum, as well as F . For AR/Rq < 1 and defining AT as the observed 
variability time (in seconds), e.g. the observed FWHM of the flare or spike in the light curve, 
we have / /AT ~ TfF e\~ a ~ F e\~ a AR/Ro, which can be used in Eq. ([23]) to determine 
L : 

Lo „ ^4 A J^L ( I±f ^ °~* £? -i = 1.3 x 10" f 4 Az« e«-i erg s" 1 , 

AT Ai2 V 2 y 1# Ai? V 2 y L ' zl AT 

(119) 

where f = 10 _6 / -6 er g cm -2 , and d,L, z i is g?£ in units of di{z = 1) ~ 2.05 x 10 28 cm (for 
standard cosmological parameters). It is hard to determine AR/Rq without a detailed fit 
to the model spectra (see Fig. [12] for the dependence of the time integrated spectrum on 
AR/Rq), and this is a price for the simplicity of this method and the use of simple analytic 
formulas rather than numerically evaluating the set of nested integrals in order to calculate 
our model predictions. One can either assume AR/R Q « 1 or try to estimate its values by 
eye, guided by Fig. [12], if one wishes to avoid a direct fit to the model predictions. 

The quantities a, L , and £i* may in turn be used to determine T"q*Rq. In order to do 
this in practice, we need to use Eq. (1551) and the relation 

(1 + z)e u = (rJCa)- 1 ^ = [249C 2 (a/2) 5 / 3 10 4 ( Q - 2 )L ^ 2 (r , 2 ) 2a J R ,i3] , (120) 
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where C a = 100C 2 is a coefficient whose value is determined numerically. The dependence 
of ei* on r , Rq, an d a is demonstrated in Fig. [13j Since a test photon of dimensionless 
energy e u pair produces primarily with photons of energy ~ T 2 /(l + z) 2 Ei*, and we fix the 
value of L [i.e. the photon number density near (1 + z)e = 1], the values of £i* becomes 
almost independent of a near To ~ a/(1 + z)eu (which is roughly where the lines for the 
three values of a for the same value of £i* almost meet). 

Eq. (I12ip can be inverted in order to obtain 



(r , 2 ) 2a iVi 3 = C 2 - 1 40.2 2 - a (a/2)- 5 / 3 L , 



r >2 



1 + z)e u m e c 2 
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If one makes the additional assumption that Rq ~ TqcAT/ (1 + z), which is valid for a large 
class of models, then Eq. (1122)1 provides the following estimate for Tq: 



T « 100 (1 + z 
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For GRBs, one may perform a consistency check for the assumption that Rq ~ TqcAT/(1 + 
z) by comparing the value of T under this assumption from opacity to pair production 
(Eq. |123j ) to the estimate for T from the time, T dcc , of the onset of the afterglow emission[] 



ro(Tdec) 



(3 - k)E iso 



l/2(4-fc) 



i ' 4B - e 'iso 1 53 n J dec,2 



131 -E iS o 53 At 



l/4 T -l/4 



dec,0 



(* = 0) 
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where Td ec = Td ec ,o s = 100Td ec> 2 s, E^o = 10 -E iso 53 erg is the isotropic equivalent kinetic 
energy in the outflow, p ext 
with radius, which is p cxt - 



AR is the external density, and is assumed to be a power law 
for a uniform external medium (k = 0) of number density 



ran, 



n 



Uq cm 3 , while A = 5 x lO 11 ^ g cm 1 for a stellar wind environment (k = 2). 
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A. Changes of variables in section 13.21 

A.l. Change of variable from (s,/i r ) to (R e ,R t ) 

Since we integrate over dQ r = d(j) r dfj, T and the integrand contains \d/j, e /d/j, r \, we can 
conveniently change variables from fi r to R e . It is straightforward to verify that, irrespective 
of the position of R t with respect to the upper bound of R e , \d\i e jd\i r \ d\x r = (dfi e /dR e ) dR e , 
with the integration over R e being performed from the smaller to the larger bound, and 
dn e /dR e > always. Furthermore, the perpendicular distance from the line of sight from 
the center of the emitting sphere to the observer at infinity, 

R± = Rt,o sin 6 t fl = R t sin 9 t , (Al) 
is constant along the trajectory of the test photon (see Fig. [2]). Thus 

d a d a j R±d8 t R±dfj. t RtdR t , . , 

s — Ha cos v t —Rt cos u t , ds = n — = 5-777 = — = ~ Rt ndR t , ( A2) 

sin 2 ^ (l-p2) 3/2 ^Rj-Rl ' 1 ^ 



7 This estimate is for the Lorentz factor of the outflow after the passage of the reverse shock, so it is close 
to that of the original outflow before it was decelerated by the reverse shock only as long as the reverse 
shock is at most mildly relativistic. For a highly relativistic reverse shock, the original Lorentz factor of the 
outflow can be much larger than this value. 
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where the last approximation holds since R± < R±, max ~ R{R±,max) /^(R±,max) = 0(R t /T) <C 
R t and the contribution from R t -C -R(-R±, ma x) is negligible for T ^ 1. 



A. 2. Integration over d(fi r and &E{ 



The local photon field derived above is symmetric around the radial direction (i.e. does 
not depend on <p r ). As a consequence, (fi r appears only in the function cr* [x( £ t, £ u l^u)] (1 — 
fiu), where fi t i is a function of cos</> r (see Eg. |21) ). Thus , we c an write f^ dcfrr = 2 J* d<fi r . 
Next, we follow the insights of (jStepney fc Guilbertlll983l ) and (jBaring)ll994j ) by performing 
the change of variables (£i,(fi r ) — » (x, u )> with x 2 — £ t £ iU defined in Eq. ( |T8l) . Defining 



C+ = [1 - cos (0 r + 9 t )\ /2 and C- = [1 
Eq. (J2U) yields /„* dfa = 2 J f C _ + du [(C+ - C 
wise Jy et dei = 2 f^xdx^tu)^ 1 - Eq. (I3"2"|) now reads: 



cos i 
sin i 



9 t )\ /2, where C+ > C- for (9 t ,9 r ) E [0,vr] 



) sin^p = 2 j£ du [(C + -«)(«- C-)] 



-1/2 



Like- 



Tjj{£t, 9tfl, Rtfi 



8<jt 



(4Tr) 2 m e c 3 R tfi 7i 



f j ~ 5 3 du 
udu 



dR e — — 
V dR P 



c_ > /(c + -u)(u-C-) A x ot 



dx°*(x) Tl 



J X 2 /etCS 



(Re) ■ 



(A3) 



A.3. L' e ,(R e ) = L' Q (e[f- a x h(R e /R ). 

i 

The specific luminosity in the co-moving frame is conveniently parameterized as L' e , = 
L' (e') 1 ~ a h(R e / Ro) where h(l) = 1 is normalized at R e /Ro = R e Rt/Ro = 1- Similarly, we 
want to parameterize the specific luminosity in the lab frame at Rq as L £ (Rq) « Loe l ~ a , even 
though the luminosity at a given radius is not really well defined, since the Doppler factor 
also depends on the angle 9t t o from the line of sight. The normalization coefficients in the lab 
frame (L ) and in the co-moving frame (L' Q ), are the specific luminosity at Rq corresponding 
to a photon energy of m e c 2 ^ 511 keV in the respective frames. Since the typical value of 
the Doppler factor is 5 = e/e' ~ T, and specifically S(Rq) ~ T(R ) = T , these coefficients 
are related by L§e x ~ a ~ T L' Q (e') 1 ~ a and L ~ TqLq. Thus motivated, we use this relation 
as the definition of L , L = T°L' . Therefore, L'(R e ) = "^(e/) 1 "" x h(R e R t /R ). It 

i 

is convenient to express the optical depth r 77 in terms of L which is approximately the 
observed isotropic equivalent luminosity at an observed photon energy of 511(1 + z)~ x keV 
near f ~ 1 for AR > R . For AR/R « f f /(l + m) < 1 (see Eq. [log] ), the peak isotropic 
equivalent luminosity at 511(1 + z)~ x keV is ~ TfL and the corresponding optical depth 
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near the peak of the spike in the light curve at the same photon energy is ~ TjTq. Therefore, 
Lq is practically an observable quantity, making it convenient to work with. 



Eq. (IA3I) now becomes: 

2r - a L erV T 



r 77 (e t , 9t t o, Rt,o) 



(A7T) 2 m e c 3 R t , 

u a du 



dR t 



dR P 



5 2+a dfi 




, x — / d X X~ (2a ' 1] ^(x) 

C- V(C + -«)(«-C-) *r ' 



(A4) 



From Eq. (8) of (IBaringl I1994J ) we can write, to a very good approximation: 



°T Jl 



7 



6a 5 / 3 



(A5) 



Making the change of variab le u — > t — {u — C-)/(C+ — C-)> we 8 ^ so have (see Eq. 15.3.1 in 
Abramowitz Stegunlll980l ): 



C+ 



u a du 



-At 



7T 



2 Fi(-a,0.5;l;-C) (A6) 



where 2 F X is a hypergeometric function and ( = (£ + — C-)/C- > 0- We define H Q (() = 
2F 1 (— a, 0.5; 1; — £), and notice that it is regularized by the factor when C_ — > 0. Note 
th at C is of order unity and H n (C) is a simple polynomial when « is an integer (see Eq. 15.4.1 



m 



Abramowitz fc Stegunlll980l . and our eqs. (1 1 1 f) - ( TT03T) for a=l,2, and 3.). 



B. Justification for the Approximations in Case 1 

When calculating the photon field at the instantaneous location of the test photon 
on case 1, where the test photons lags behind the shell, R t < R s h(t t ), we have used an 
approximation for the value of 8 r , namely Eq. floT?]) . which is valid for 9 r 1 and break 
down for 9 r ~ 1 which corresponds to 1 — R e — 0{Y^ 2 ). This is despite the fact that in 
this case 9 r can assume any value between zero and it. The justification for this convenient 
approximation is that the contributions to the optical depth r 77 from 9 r ~ 1, where our 
approximation breaks, is negligible compared to the contribution from 9 r <ti 1, where our 
approximation is valid. In order to show this more explicitly, we examine the dependence of 
the integrand in the integration over dR e on the value of 9 r in the range 

_ < 9 r < 1 ^ < R sh - ^ e w 1 - R c < 1 , (Bl) 
1 1 1 1 
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which gives us a handle (up to factors of order unity) on its dependence throughout the 
entire range of possible 9 r values. In this intermediate range of 9 r values, i? s h — R e ~ 1 — R c 
since for case 1 T 2 (R sh — 1) < a few, and thus 



f 2 



Rsh — Re 4 



2 \ R s h — R e ) ( - 



(m + i)r 2 



+0 (I7 4 ) » (i? sh -i? c ) 2 » (l-i? e ) 2 . (B2) 



Likewise, Eq. fl50l) yields 
djjb e 1 



di? e 2r 2 i?2 



fm — R t 



r? (X\ - i ) 

m+l 



l+mR" l+l ) -R™ +2 
m + 1 v 1 



fm-l 



2{m+l)T 2 t Rl 2(m + l)r 



2 ' 



(B3) 
(B4) 



Thus we see that in this intermediate range of 9 r values d\x e jdR e is approximately constant, 
and is of order T^" 2 . Besides, H a is of order unity, ~ 9 2 , and 5 ~ T~ l 6~ 2 , so that 



S a+2 dn 



f 2 dR, 



1 



oc 



1-iL 



(B5) 



Moreover, Eq. (|69l) . which is valid for case 1 in the limit 9 r <C 1 that applies in our inter- 
mediate regime, implies that (1 — R e )9 2 is approximately constant in this range of 9 r values 
given by Eq. (1B1 j) . Therefore, from Eq. (1B5I) we conclude that the integrand in the inte- 
gration over dR e is approximately constant over this range in R e , which is of interest here. 
Furthermore, the integrand must still have a similar value, up to a factor of order unity, 
even for 9 r ~ 1, since the approximation of 9 r 1 breaks only marginally, rather than very 
severely (since 9 r cannot have values ^> 1). As the region where our approximation breaks, 
9 r ~ 1, corresponds to 1 — R e — (D(T^ 2 ), i.e. a range of the order of T^ 2 in R e , which is 
much smaller than the range over which our approximation is valid, and is also much smaller 
than the range in Eq. (IBlj) . we conclude that the contribution to the integral from 9 r ~ 1 
can safely be neglected. 



C. Properties of the Photon Field in Case 3 

By changing the integration variable from fi r to R e we eliminated the need to express 
R e as a function of fi r , and to calculate the minimal value v r which corresponds to # rimax - 
Nevertheless, this is still interesting in terms of the properties of the local photon field, so 
it is given in this appendix. Each value of [i r may correspond to two different values of R e , 
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one at the front and one at the back of the equal arrival time surface of photons to the point 
(Rt,t t ). Eq. fl86l) can be re-written as 



R. 



m+2 



iC+i + 2(m + l)V 2 t (l - fi r ) R e + 2(m + l)r 2 (l - // 

'ct 

For m = this becomes a second order equation with the solutions 



e,max 

R™ +2 - 2(m + l)r 2 ( ^ - ^ r ) i? e + 2(m + l)r 2 (l - // r ) = . 



(CI) 



i± Wi 



2r 2 (i-/y) 
r^K/i^-i; 



i ± 



\ 



i - 



4(r t r 



-Re.max + (Tt0 r ) 



(C2) 



where ^ max may be obtained by the condition of a single solution, i.e. that the expression 
in the square root vanishes. This implies 



2-R P 



\ 



Rp 



2 — R P 



(C3) 



where we chose the root of the equation which corresponds to the familiar result of T t 9 rjUiax ~ 

Re, max/2 for Re : ma,x !■ 

More generally, fi r ,mm = cosd r)Uiax may be found by the condition that dfi r /dR e = 0. 
Using Eq. (IBB]) this results in 

1 m+2 



m + 1)^ +2 - (m + 2)iC +1 + R m+1 







e Vr.max I 



(C4) 



Alternatively, one can use the latter relation, which is obtained by substituting dfi r /dR e = 
from Eq. flBB|) into Eq. flB6l) . to obtain an explicit equation for ^ max , 



{m + l){T t 9, 



(m + 2) (r t r , r 



>2(m+l)/(m+2) . pm+1 



_i_ pm+i _ q 



(C5) 



D. The Scaling of r 77 with T 

It is instructive to explicitly derive the scaling of r 77 = To(et, Rt)J-(x) with T, in the 
three regimes 1 ^> T < Tf, 1 <§C T < Tf, and T ^> Tf. The only time dependence of To on T 
is through R = y~\T /Tq)- 1 /^^ (see Eq. [38]), so that r oc (1 + f)^ 1+am ^^ m+1 \ 
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D.l. 1 > T < T f 

For 1 ^> T < Tf, r is thus approximately constant and the time dependence of r 77 is 
dominated by the the time dependence of J-"(x), which we now consider in more detail. First, 
the maximal value of the emission angle 9 tt0 and correspondingly of x — (7^0^,0) 2 along the 
EATS-I is given by 



T-Tn 



R0O 2 



t,o 



2c 



T (m + l)x t 



•J'max ('^') 



T 



(m + 1) 



(Dl) 



This result holds in general, and can be readily obtained by noticing the x max always cor- 
responds to y m - m = (T '/To)~ 1 / ( - ,n+1 \ and substituting the latter in Eq. (TI5]) . Therefore, 
x < x max <C 1 for T <C 1 and we are always in case 3. Second, it is straightforward to show 
that 



fm(Rt 



R + AR 



-(m+l) 



n, ) \i + T f 

so that for T < Tf we always have 



1 + T\ 1 + x(m + - Rt l ) 1+T 



1 + x(m + 1) 



and 



i?e, max = Re,3 = fm(Rt) 1/(m+1) = R^ k + x(m + 1) ( 1 - i^ 1 



1 + T* 



l/(m+l) 



(D2) 



/■-Re.3 

<iR t / dR e T(R t: R e ,x) 

1 J Ro/Rt 



(D3) 



(D4) 



where X(-Rt, i? e , is given in Eq. ( 1971) . and 



6 _ _ d-1 

-'•'e.niin * 



1 + (m + l)x 
1 + f 



l/(m+l) 



(D5) 



Keeping terms to first order in T (and x), the range of R e value that is being integrated over 
in Eq. (|D4]) is 



AR P = R, 



e.max ± ^e,rmn 



Re 



T 



x 



{m + l)R t R: 



^ = 0{T) « 1 . 



(D6) 



The integrand includes in several places the expression 



f m (Rt) - RT 1 = R m+l 



pm 

pm+l ^ e.max , R 
e,max ~ /.__ , i \ \ ^ejnax 



(m + 1) 



Re. < 1 



(D7) 
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which is either comparable to or much smaller than 1 — R e , which also appears in the 
integrand, thus defining different regimes. The relevant ratio to compare to unity is 

t-Re,max — -Re \ AR e 1 — -Re.min 1 / A i\~ 1 f T x\ . . 

— ' ~ = ~ = — 1 ~ [Rt — 1 7 — — , (D8) 

I -Re J 1-ReM* V J \ m + 1 Rt) 

which measures both the fractional change in 1 — R e and the minimal value of its ratio to 

-Re, max -Re- 

For x = x max = T/(m + 1) this ratio is x max / ^t<f/(m+l)<l so that 1 — R c is both 
approximately constant and much larger than -R ei max — -R e ~ fm(Rt) — -R™ +1 . Therefore, the 
only term that varies significantly with R e in the inner integrand is (^H a (Q. For £ ^> 1, 
H a {C) ~ C" so that CH a (Q ~ (C-C)° oc f a / 2 (R £jmax - R e ) a ' 2 where the integration over 
(-Re,max — R e ) a ^ 2 results in a factor of T 1+a / 2 so that altogether the inner integral is oc T a+1 . 
The outer integral is of the form J" 00 dR t g(Rt) = const. For £ < 1, H a (Q ~ 1 and (°H a (() ~ 
~ (r^ r — r^ 4 ) 2a which consists of a sum of terms of the form T a (i? ejmax — R e ) a ~ a that 
upon integration are oc T a+1 . Thus, 

^(x max ) oc f Q+1 . (D9) 

Note that in this case most of the contribution to the optical depth comes from R t < 2 or 
AR t ~ 1. 

For x = 0, 9 t = ( = so that H a (() = 1. Furthermore, the ratio in Eq. (lD8h becomes 
larger than unity for R\ — 1 < T/(m + 1), and in this regime f m {Rt) — R™ +1 ~ -R e ,max — -R c ~ 
1 — _R e so that is roughly constant and the inner integrand scales as (1 — -R c ) _1 , which 
upon integration scales linearly with T, 



/ dRt g(Rt) / — « 0(1) / d(& - 1) In 



'o--*msx.y&t (i - -R 



(At - 1) + x B 
(Rt - 1) 



<7(l)(21n2)x max ocT. (D10) 



For R t — 1 ^> x max = T/(m + 1), the approximation discussed in the previous paragraph 
apply, and this part of the integration over R t does not contribute significantly to the total 
optical depth, so that 

J r (s = 0)ocf. (Dll) 

Physically, the lack of significant contribution to the optical depth from R t — 1 ^> x max = 
T j (m + 1) may be understood since the maximal value of 9 r (which corresponds to R e = Rq) 
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starts to decrease significantly, 
max[(r # r ) 2 ] 



(m + 1) 



R t 



Rt-l + 



T 



n -1 



(m + 1) 
4-K 



(m+l)(H t -l)flf 

which suppresses the opacity to pair production. 



(m+l) ' 
«1 ^-l»(^ T y, 



(D12) 



D.2. 1< T < T/ 

For 1 <C T < Ty, we have i? i ^ Ro, so that Ro/Rt *C 1 and may effectively be taken 
as zero. Furthermore, -R eima x < Rl{T) < R L (Tf) = R + AR (since X < Xy) so that i? e>2 = 1 
and R e>3 = f m (x, R t ) 1 ^ m+1 ^ is given by Eq. (ID3I) . and Eq. ([9"oD now reads 

fR 2 {x) pi _ roc rfm(x,R t ) 1/(m+1) 

T{x) = dR t dR e I(Rt, R e ) + dRt dR e T(Rt, R e ) . (D13) 

Ji Jo Jr 2 (x) Jo 

In this regime neither the boundaries of integration nor the integrand, X, depend on T. As 
a consequence, the dependence of r 77 in this regime is only through r , and we have 

r 77 (l < f < f f ) w r (f ) J-(x) oc f N-^fl/M _ (D14) 



E. On the definition of the optical depth 



We start with the explicitly Lorentz invariant expression for the differential interaction 
rate of two particles, denoted '1' a nd '2', colliding with respective momenta pi and p 2 , as 



given in Eq. (24a) of I Weaver! (119761 ) : 



2™ 2^41 1/2 



Ru(Pi,P2) 



ni(pi)n 2 (p 2 )(l - Pi ■ 2 ) [(Pi -P2) ~ 

Pi -P2 



-ca 



(El) 



where pi, p 2 are the four-momenta of particles '1' and '2', respectively, mi and m 2 are 
their masses, ni(pi) and n 2 {p 2 ) their phase-space density and a is the generalized Lorentz- 
invariant cross-section, usually computed in the center of momentum frame. In Eq. flEll). we 
have e xplicitly written the dependence of R\ 2 on the momenta, which is missing in IWeaver 
( 119761 ). in order to distinguish it with the total interaction rate < -R12 >• The latter results 
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from an integration over the phase spaces of both particles (see eqs. (2), (27) in ([Weaver 
1976h ) : 

< Rn >= Y^~ 2 J J R u(PhP2)d 3 pid 3 p 2 . (E2) 

In Eq. (|E2|) . the Kronecker symbol 8x2 is 1 if the two particles are identical and otherwise. 
It accounts for the fact that, for identical particles, the double intergration counts twice each 
pair of interacting particles. 

Now, we define R\2{p\) as the interaction rate of a given particle '1' of momentum p\. It 
writes Rnipi) — f Ri2(Pi,P2)d 3 p2, without a Kronecker symbol because there cannot be any 
double counting when there is no double integration. Specializing now to 77-interactions, the 
interaction rate i? 77 (pi) is equal to the decrease in ri\ per unit time : drii(pi)/dt = — i? 77 (pi). 
Defining the optical depth of a particle of type 1 and momentum p\ as the corresponding 
attenuation per unit length : dni/rii = —r(pi)ds, where ds is an element of trajectory of 
particle 1, we obtain : 



T 77 (pi) = -R 77 (pi)/cni = / n 2 (p 2 )(l - Pi ■ P2)crd 3 p 2 



(E3) 



where in the last equality we made use of mi = m 2 = in Eq. flEll) . We thus re-derived 
Eq. ffl6|) (in integral form), and showed that there is no factor 1/2 involved because the 
computation of the optical depth does not warrant a double integration over the phase space 
of both partic les. Because they comput e the t otal reaction rates and not the optical depth, 
Weaver! (119761 ) and IStepney fc Guilbertl (119831 ) do have this factor. 



Another source of confusion arises from the fact that in their seminal paper. iGould fc Schreder 



(119671 ) specialized to an isotropic distribution for particles 2, which brings up a factor 1/2 
due only to the normalization of the integration over cos#. In other words, introducing 
dn = n 2 (p2)d 3 p2 = (l/2)n(e)de sin 6d0 in Eq. (lE3h immediately yields their Eq. (7). 
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1 12 * j_ * 
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Eq./§ 




ooseiveu pnoion eneigy iioimanzeQ uy i iic election iest eneigy 
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spherical coordinates (time, radius from the source, polar angle) 




i?o, AR 


onset radius and range of the emission episode 
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radius and bulk Lorentz factor of the emitting shell 
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— dlogR 
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power law index ot co-moving spectral luminosity with radius 
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test photon spherical coordinates and shell Lorentz factor at Rt 


§S 
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initial test photon spherical coordinates and Lorentz factor 


mm 




emission radius and polar angle of interacting photon 




T ,T 


arrival times of first and subsequent photons at the observer 
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maximal radius ot emission along the LAI b-1 and LAI 
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cosinp of angle 8 
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x = y 2 


i*x*i ii *j_i j_r j_r 

dimentionless photon energy m the center ot momentum frame 




C = (1 - A*«)/2 


convenient integration variable 
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r 


interacting photon emission to test photon intersection distance 


§13 


r 


angle of an interacting photon relative to the radial direction 


§E1 


5 = (l + z)e/e' 


Doppler factor between the co-moving and lab frames 


m 


fm 


useful quantity 


(J62D 




typical optical depth at e = 1 on a dynamical time (Tf > T ~ 1) 


(J35D 


To, Fix) 


explicit analytic and integral parts of the optical depth 


(EHD 


x = (r t , ^,o) 2 


rescaled emission angle squared 


(USD 


y = Rt >0 / R L {T) 


1* 1 1 i xl • 1* i — \ Amn T 

emission radius rescaled to the maximum radius on an LAlb-1 


/loll 


R = Rj Rt 


radius rescaled to a given test photon emission radius 
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R = i?/ i?j 


radius rescaled to the instantaneous test photon radius 
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f = T/T - 1 


arrival time of photons rescaled to the earliest arrival time 


§0 


5 = S/T t ,fl e = T 2 t fx e 


rescaled Lorentz factor and cosine of the emission angle 


(EHD 


C- = rfC- 


rescaled integration variable 


(EED 


.(/max //mm 


rescaled variable y, and value at which T changes its behavior 


(imp 



Table 1: Notation and definition of some quantities used throughout this work. 
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m 
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log 10 (A J R/ J R ) 


Figures 


2 








-2,-1,0,1,2 




2 





1 


-2,-1,0,1,2 


m 


2 


3 


-2 


-2,-1,0,1,2,4 


\m 


2 


0,1,2 





2 


M 


2 





-2,-1,0 


2 


M 


2,3,4 














Table 2: The different sets of parameters for which results are shown in this work. 
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test photon case A case 2 case 3 




Fig. 1. — An illustration of the two different equal arrival time surfaces (EATS) of photons: 
1. to the observer at infinity (EATS-I, in blue), and 2. to the instantaneous location of a 
test photon (EATS-II, in different colors). The overall geometry as well as relevant radii and 
angles are shown in the upper panel, along with an illustration of the three different cases 
that are discussed in the text, in which the test photon either lags behind the shell (case 
1), coincides with the shell (case 2), or is in front of the shell (case 3). The lower panel 
shows the sequence of EATS-II, whose size increases with time, nested within the EATS-II 
which correspond to a larger time, and in particular within EATS-I which corresponds to 
an infinite time (i.e. an infinite radius for the test photon, when it reaches the observer at 
infinity) . 
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Fig. 2. — Geometry of the interaction between two photons, for a spherically symmetric 
shell. A test photon emitted at Rt$ reaches Rt > Rt,o at t t > to and can interact with a 
photon emitted at R e that reaches the location R t at the exact same time tt as the test 
photon. Note that O, Rt$ and Rt are coplanar (and in the plane of the figure), whereas R e 
is not in the same plane, nor is the interacting photon trajectory that goes from R e to R t . 
The observer is to the right, at infinity. The other symbols are defined in the text. 




Fig. 3. — The equal arrival time surface (EATS-II) of photons to (Rt,t t ), which represents 
a general point along the trajectory of a test photon, is shown by the thick blue line. It 
naturally divides into three cases: 1. the test photon is behind the shell (R t < R s h(t t ) - 
upper panel), 2. the test photon coincides with the shell (R t = R s h(tt) - middle panel), and 
3. the test photon is ahead of the shell (R t > R s h{tt) - lower panel). There are qualitative 
difference in the properties of the EATS-II between these different cases, that are discussed 
in the text. 
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Fig. 4. — f m (R t ), defined in Eq. fl62l) . as a function of R t , for m = and m = 1. The test 
photon is necessarily on the shell at the time of its emission, so that all the curves meet at 

/m(l) = 1- 
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10 10 10 10 10 



x 1/2 = YtAo 

Fig. 5. — The normalized optical depth F(x) = r 77 /ro, as a function of the renormalized 
emission angle, x 1//2 = 'jt.odtfl, for several different emission radii -R^o- The upper panel is for 
AR/Rq = 1000 while the lower panel shows the results for AR/Rq = 1 (solid lines) and for 
AR/R = 1000 (dashed lines) overlaid on each other. The small vertical lines in the lower 
panel indicate the angle that corresponds to T = Tf, outside of which the contributions to the 
opacity from R > R + AR for AR/Rq = 1 start being missed (this effect becomes significant 
only at somewhat larger angles; see discussion in the text). For i? ti0 = 2R = R + AR, this 
corresponds to x 1//2 = 7^0^,0 = 0, which is outside the range shown in the figure. In both 
panels the photon index is a = 2 while the Lorentz factor and the total luminosity in the 
comoving frame are independent of radius (m = b = 0). 
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Fig. 6. — The normalized optical depth J-{x) = r 77 /ro, as a function of the renormalized 
emission radius, (Rt t o/Ro) — 1, for different values of the normalized emission angle x 1 ! 2 = 
7t,(A,o- The upper panel is for AR/Rq = 1000. The lower panel is for AR/Rq = 1 but also 
shows the corresponding result for AR/Rq = 1000 in dashed lines, where the x-symbols show 
the value of the emission radius corresponding to an observed time of T = Tf (for 7^0^,0 — 1 
this corresponds to (Rt,o/Ro) — 1 = which is outside the range shown in the figure). Note 
the deviation near (R t fi/Ro) — 1 ~ 1 and see the text for discussion of its origin. 
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Fig. 7. — The normalized optical depth = r 77 /r , along the equal arrival time surface 
of photons to the observer (EATS-I), for several different values of the normalized time 
T = (T/T ) — 1: the upper panel shows J-{x) as a function of the normalized emission radius 
y = R t fll Rl(T) for several values of 1 > T < Tf, while the lower panel shows ^{x) as a 



function of Y = (y - y min )/(y n 



[ 0C r 



i)/i max for several values of 1 <C T < Tf. 
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Fig. 8. — Lightcurves (left panels), instantaneous (thin lines) and time integrated (thick 
line) spectra (right panels), calculated using our semi- analytic model, for a constant Lorentz 
factor (m = 0) and a comoving emissivity independent of radius (b = 0) with equal energy 
per decade of photon energy (corresponding to a photon index of a = 2). We show results 
for three different radial extents of the emission region, AR/Rq = 0.01, 1, and 100, from top 
to bottom. We also use r* = 1 (see Eq. [35]). 
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T = T/T -1 e [m e c 2 ] 

Fig. 9.— Similar to Fig. El with b = 0, a = 2, r* = 1, but for a fixed AR/R = 100 and 
varying m where T 2 cx R~ m . 
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T = T/T -1 e [m e c 2 ] 

Fig. 10.— Similar to Fig. El with m = 0, a = 2, r* = 1, but for a fixed AR/Rq = 100 
and varying b where the the spectral luminosity in the comoving frame of the shell scales as 
L' s , oc R\e') l - a . 
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Fig. 11. — Similar to Fig. [HJ with b — m — 0, t* = 1, AR/R = 1 and varying a, where the 
the spectral luminosity in the comoving frame of the shell scales as L' £ , oc R h (e') l ~ a . The 
middle panel and bottom panel are for a = 2 and 3, respectively. The top panel is for a = 1, 
for which r 77 becomes independent of the photon energy e and therefore the spectrum is 
always a pure power law, F e oc e° and the flux depends only on time but not on the photon 
energy. For this reason we show light curve (left panel) and time integrated spectra (right 
panel) for different values of r* (see Eq. |35j). 
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Fig. 12. — The time integrated spectra for different values of AR/Rq, while fixing the values 
of the other parameters. In the top panel m = b = 0, in the middle panel m = 3 and b = —2, 
while in the bottom panel m = and 6=1. The values of the other model parameters in all 
the panels are a = 2 and r* = 1. 
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contours of the spectral break energy £-|*m e c 2 
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Fig. 13. — Contour lines for of the photon energy 6um e c 2 where the time integrated spectrum 
over a flare or spike in the light curve steepens due to opacity to pair production, shown in 
the T — Ro plane, according to Eq. (j!2ip . In the upper panel ei 1f m e c 2 is varied and a = 2 
is fixed, while the lower panel a is varied and Eum e c 2 = 1 GeV is fixed. In both panels 
-^0,51 — C*2 — z — 1- 



